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ABSTRACT 


As the accuracy of the technique of long-baseline radio interferometry 
improves, the errors contributed by the extended structure of natural sources 
will become increasingly important, fo begin a study of such structure effects, 
this report presents a theoretical framework, proposes an effective-position 
approach to structure corrections based on brightness distribution measurements, 
and analyzes examples of annlytical and measured brightness distributions. Other 
topics include the effect of the frequency dependence of a brightness distribution 
on bandwidth synthesis (BWS) delay, the determination of the absolute location of 
a measured brightness distributiovi, and structure effects in dual-frequency cali- 
bration of charged-particle delays. For the 10 measured distributions analyzed 
in this report, it is found that the structure effect in BWS delay at X-band (3.0 
cm) can reach 30 cm but typically falls in the range of 0-5 cm. The largest BWS 
delay effects occur at points in the u-v plane with very small fringe visibility. 
To provide a possible method for reducing and estimating structure effects in BWS 
delay, a limit approach is investigated through a general analysis and through a 
study of particular analytical sources. A trial limit equation that is dependent 
on visibility is successfully tested against the 10 measured brightness distribu- 
tions (seven sources) . If the validity of this particular equation for an upper 
limit can be established for nearly all sources, the structure effect in BWS 
delay could be greatly reduced without supplementary measurements of brightness 
distributions, perhaps to the 2 cm level (lo) on intercontinental baselines and 
to less than a centimeter on regional baselines. 



I. INTRODUCTION 


As the accuracy of the technique of long-baseline radio interferometry 
improves, the effects of source structure will become more Important in measure- 
ments of baseline, source positions, and other quantities. A crude method for 
demonstrating the magnitude of structure effects is to note that changes in source 
position due to structure are expected to be' of the order of the fringe spacing, 
which is given in radians by the ratio X/B where A is the observing wavelength 
and B is the baseline length. Since the partial of delay with respect to source 
position is of the order of B in meters /radian, the corresponding effect on delay 
will be of the order of (A/B)B s' A, Thus, for observations at X-band (A ■ 3.6 cm) 
some sources will contribute unacceptable structure errors if centimeter-level 
accuracy is desired. 

In addition to the geophysical applications of VLBI, differential VLBl mea- 
surejnents con be used to measure the separation between close pairs of radio 
sources (natural-natural and natural-spacecraft). Source structure eff»»f’ts will 
he even more significant in this application since error sources other than 
source structure are expected to contribute a total error of only 10"^ to 10"^ 
arcseconds in position on intercontinental baselines. This potential accuracy 
is to be compared with the estimated structure error (A/B), which is equal to about 
10“3 arcseconds at X-band with B “ 6000 km. Thus for differential measurements to 
reach their full potential, source structure effects must be treated very care- 
fully. 

To begin the investigation of source structure effects, this report presents 
a theoretical framev^ork, proposes an effective-position approach to structure 
corrections based on brightness distribution measurements, analyzes a few examples 
of analytical and measured brightness distributions, and develops a limit formula 
for reducing and estimating structure effects in bandwidth synthesis (BWS) delay. 

The report is organized in the following way. The next section, Section II, 
rederives the cross-correlation function by beginning with a description of the 
noise wave emitted by a natural source, and tracing the signal through the 
instrumentation and cross-correlation process. The derivation is repeated in 
detail to moke this report self-contained and to emphasize the relation of the 
parameters of the brightness distribution to the parameters of the cross-correlation 
function. That relation is needed to demonstrate clearly how the absolute location 
of an independently measured brightness distribution can be determined through 
standard VLBI measurements. Section III reformulates the brightness transform to 
cast it in a form mote suitable for the analysis of the effect of structure on 
VLBI observables. Section IV discusses the determination of structure effects in 
VLBI observables (BWS delay, delay rate, phase delay) and introduces the concept 
of effective position. Sections V and VI consider the structure effects associated 
vi^ith bandwidth synthesis (BWS) delay and with delay rate, respectively. Section 
VII combines the results of Sections V and VI to derive the effective position 
associated with BWS delay and delay rate. Section VIII considers the structure 
effects associated with the phase-delay observable. Section IX shows how to 
correct delay and delay rate for structure effects and in the process demonstrates 
a technique for determining the absolute location for the origin of the coordinates 
with which the brightness distribution is specified. Section X considers differ- 
ential measurements of the relative positions of two natural sources and notes 
that such measurements might be used as a direct and accurate test of the proce- 
dures used to correct structure effects. The subsequent four sections present 
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examples of brightness distributions to illustrate the general analysis. Speci- 
fically, Sections XI, XII and XIII treat the analytical e.tamples of a symmetric 
distribution, a double-point source and an N-polnt source (with emphasis on triple 
point), respectively, while Section XIV presents results for actual natural 
sources. Finally, Section XV derives and tests a limit formula for structure 
effects in BWS delay. 


II. THE CROSS-CORRELATION FUNCTION 


This section presents a derivation of the cross-correlation function for an 
extended natural source that is completely incoherent. The derivation begins 
with a description of the noise wave emitted by the source and traces the signal 
through the Instrumentation and cross -cor relation process. Considerable detail 
is Included so that the parameters of the brightness distribution are precisely 
defined in relation to the cross-correlation function, Such care is required to 
show chat it is possible to use the standard multiparameter fits to delay to 
cetermine the absolute location of the unknown origin of an independently 
measured brightness distribution. As the following derivation shows, a distinc- 
tion must be made between the actual effective position, Che origin of structure 
coordinates and an assigned reference position (see Section IV). 


The radio noise generated by a very distant extended natural source can be 
expressed as a superposition of plane waves integrated over possible wave direc- 
tions and frequencies: 


E(x,t) 


CO ^ 

If + c.c. 

k 0 


( 1 ) 


where E(x, t) is the electric field at point x and time t,A(k,oj) is the random 
Fourier amplitude at frequency w for the wave portion received from direction 
(c, and the letters c.c. denote complex conjugate. The wave is assumed to be 
linearly polarized for simplicity. All quantities are measured with respect to 
a quasi-inertial geocentric frame with axes defined by true equatorial coordinates 
of date. The wave direction k, which includes annual aberration, must be ex- 
pressed as a function of two parameters. In terms of right ascension and de- 
clination, the wave direction is given 

A 

Ic » -(cos 6 cos a, cos <S sin a, sin 6) (2) 


where a, 5 are apparent right ascension and declination relative to true 
equatorial coordinates of date. The quantity dfi represents a differential 
solid angle for integration over the possible wave directions and is^glven by 
cos 6dad6 in the case of right ascension and declination. The term k in the 
argument of the Fourier amplitude stands for the two direction parameters. 
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The eleccjric field detected at antenna j is given by 


Ej(t) - E(x^(t)»t) 


( 3 ) 



k 0 


A(k,w)c 


iwft-k'x^ (t)/c] 


dwdj? + c.c. 


(A) 


where iL (t) is the location of antenna j at time t. The voltage signal recorded 
at antenna j Is given by the expression (Ref. 1) 
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( 5 ) 


where G. is the effective system bandpass filter and n. is random instrumental 
noise. ^The phase represents all Instrumental p^vise^ effects and includes the 
heterodyne phase, cable delays and bandpass phase shifts. For simplicity, we 
will neglect transmission media effects since they will only add more terms to 
the phase and will be of no significance in the derivation of the cross-correlation 
function. We have assumed that the antenna pattern is large compared to the 
source size and can therefore be neglected. 

In the following, we will neglect two^level sampling and treat only the analog 
signals. As is well-known, the digital result la the same as the analog result 
in the case of weak signals, except for a loss of 2/tt in SNR. Further, we will 
neglect clock synchronization errors and will not treat delay offsetting of both 
signals. These considerations modify the phase only slightly, can be readily 
analyzed and corrected with adequate a priori Information, and have no ultimate 
effect on the brightness transform derivation. 


In the cross-correlation procedure, the signals from the two antennas are 
offset by a model delay and multiplied together, which gives 
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where 




^ “ (w^w')t 4 - w*k’*Xj(U + %)/c - wk*x^(C)/c 4 *'}<^(t H- ^,^) “ %Ct) ■• ( 7 ) 

and (}>^ is a sdmilar expression that will not he needed. To obtain the mean 
vain 0 “ for the voltage product i one can compute the ensemble overage; 


<V^(t)V^(t 4- y> - 
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If 1 1 <A(k,m)A(k’,u5’)> e^'^cdwdfidw'dS^' 


k' 0 k 0 
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Since the instrumental noise i?j uncorrelated between antennas and Is nncor-^ 
related with the radio noise, all terms involving the additive Instrumental noise 
have averaged to zero. We will assume that the source ie completely incohey’ent 
so that the source signal from one part of the source (k)^is completely uncor- 
related with the signal from another part of the source (k’). Further, we will 
assume that a signal component emitted at one frequency (w) is completely un- 
correlated with a signal component emitted at another frequency (w'). These 
assumptions can he represented mathematically by a condition on the ensemble 
average of the component product: 

<A(k,u))A*(kIw’ )> ■* i DCk,w)6(k - k')6(w - w’) (9) 

where D is the brightness distribution at frequency m and 6 (z) represents a 
Dirac delta function. (The delta function 6(k - k') must be parameterized in a 
manner consistent with the differential solid angle dQ, but this detail will be 
transparent to us). 


h 


Since the signal is teal, we know 


* w * 

A(k,u) ■ A (k> - u) 


( 10 ) 


feo that 


<A(k,w) A(k',w’)> 

- |D(k,w) d(k - k') 6(u + w’) 

(11) 


* 0 for w, la’ > 0 

(12) 


Wnsn theae conditions are substituted in Eq. (8) we obtain 


<V^(t)Vj<t + T„)> 


/ 


A i^ -ip 

D(k,w) G^Gj e ^ doidH + c. c. (13) 


k 0 


where the phase is given by 


^ »V 

<{ig *• wk'B^/c - l|/^(t) + !|/^(t + tj^) “ WTjjj (l/») 

for which the retarded baseline is defined by 


» Xj (t + tj^) - x^(t) (15) 

The second term in Eq. (8) has dropped out since w ■ -w' is not covered in the 
region of Integration. 


To cast this result in more familiar form, we must define special direction 
parameters for Ic. The most common approach is to define two new coordinates to 
replace right ascension and declination (a, 5), where the new coordinates are a 
linear representation of small angular displacements relative to some point near 
or within the source. These coordinates will be called structure coordinates and 
are defined by the expressions 


8 ~ (c-a^) cos 6 ^ (16) 

Y s 6 - 5^ (17) 
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where (ct i 6 ) Is the location of the origin of structure coordinates in terms 
of right ascension and declination* Usually when VLBI measurements of structure 
are made, it is necessary to make strnie relatively arbitrary assignment of the 
origin of structure coordinates (e.g,, the center of the brightest component), 
Further, the absolute location of this origin is not accurately known in (a, £) 
coordinates. For the moment we will proceed as though this origin uncertainty 
did not exist and will leave discussion of its determination for Section IX. 

With the definition of structure coordinates in Eqs. <16) and (17), the differ- 
ential solid angle becomes 


dn *■ cos da di5 

r 

(18) 

« dgdY 

(19) 


provided the displacements (g, y) from the origin are small. 


Up to this point, no precise definition of source location has been made, 
since the preceding definition only established a local coordinate frame. In 
Eq. (13), the positional effects of the source enter the expression as an inte- 
gral over wave direction and no particular direction is the source "position,'' 
To proceed with the analysis, it is useful to pick some reference direction so 
that a small angle approximation can be made for the wave direction ic. Suppose 
some adopted point (3 q, Yq) within the source has been chosen as a reference 
point with a reference direction given by 


^0 “ 


If the distribution is sufficiently narrow, we can approximate the wave direction 
by 


r a. 9k 
^0 8g 


(g - gQ> + 


(Y - Yq) 


( 21 ) 


where the partlals are evaluated at (g^, Yq) * It ia important to note at this 
point that the assignment of a particular point in the source as a reference 
point is quite arbitrary. Further, it should be emphasized that the point 
(gOiYo) ia uot necessarily the same point as the origin of structure coordinates 
CSriYt) as is often implicitly assumed. Using the last relation, we obtain for 
the cross-correlation function; 


<V,(t)Vj(t+T^)> 




00 



R(u,v,w)Gj^Gje 


du + c.c. 


( 22 ) 
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where the phase is given by 


♦a ■ “''O • + T„) - (23) 

In this expression, the brightness transform is defined as 


R(u,v,o)) 



00 



D(3jY»w) 


-27ri[u(0-3«)+v(Y-Yo)3 

e u u 


(24) 


where 


A 



(25) 


A 



(26) 


Note that the limits of the (6, y) integration have been extended to infinity 
under the assumption that the narrowness of the brightness distribution 
terminates the integration. 


Using the standard definition of the geometric delay, we can write 



where Tg is the geometric delay for the direction and where v, is the velocity 
of station j . The phase of the cross-correlation function becomes 


<}>a = w(tg - T^) - if/^(t) + \p^(t + Tj^) ( 28 ) 

In the standard VLBI measurements, the first term mxg Is the quantity of Interest 
and the others are removed in processing (by phase calibration, etc.) to the 
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extent allowed by the accuracy of the calibration techniques and a priori know- 
ledge. (Had we been more thorough and considered clock synchronization, an 
additional unknown, t , would have been added to Tg. Such a term is present in 
actual VLSI experiments and its measurement is of Interest to tlme-and-frequency 
specialists.) After such corrections, the measurement of delay is still corrupted 
by the implicit phase of the complex brightness transform R in Eq. (22). It is 
the purpose of this report to Investigate this phase effect of structure. 


III. REFORMULATION OF THE BRIGHTNESS TRANSFORM 


This section reformulates the standard brightness transform in Eq. (24) to 
place it in a form more suitable for subsequent analysis. We will neglect the 
fact that the fringes are actually Integrated over a narrow bandwidth (see Eq. 22) 
and will place the observing frequency w at the center of the channel passband 
(2 MHz for the BLKO system). This will be a good approximation as long as the 
brightness distribution changes very slowly as a function of frequency over the 
channel passband. Note that a frequency change of 2 MHz would be only a change 
of one part in a thousand at S-band (2.3 GHz). Given the earlier definitions in 
Eqs. (16) and (17), one can easily show that the vectors -9(1/33 and -Dfc/Dy are 
orthogonal unit vectors in the direction of increasluti right ascension and 
declination, respectively, and foj-m a set of basis vectors^in the plane of the 
sky at the source. We will^call these unit vectors i and j, respectively, and 
will make them the basis vectors for the (3, y) coordinates. Note that the basis 
vectors at the origin (a^, 6 ) used to define (3, y) coordinates in Eqs. (16) and 
(17) will not necessarily be^exactly the same as those at (3^, Y_) • However, we 
will assume that these two points are very close and the slight differences in 
basis vectors are negligible. The brightness transform now becomes 


R = e 




If 


-iU3? *?/c 
t)e ® dn 


where dfl is the area differential and 


(29) 


•-w ^ 

P = 31 + Yj 
B = B-i + B j 

s 3 Y 


A time variable has been added to the argument of the brightness distribution to 
allow for time variability. The vector 1 is obtained by projecting iS" onto the 
plane of the sky. (Though it is not necessary to do so, we hS'/e dropped the out- 
of-plane component of to emphasize that the only in-plane component is important 
as a result of the dot product with ? - We have extracted the multiplicative 

phasor from the integral since it does not depend on the integration variables. 


(30) 

(31) 

(32) 


8 


The phase chanige (structure phase) contributed to fringe phase by the bright- 
ness transform will be given by 




(33) 


where 


e^D(?,.,t) cog (w^g* ’^/c)dfJ 
=jj D(f,w,t) sin(wig* ?/c)dn 


(34) 


(35) 


The integer-cycle ambiguity problem encountered in the computation of structure- 
phase involves some subleties that should be mentioned. As discussed in the 
next section, the observables of BWS delay and delay rate are approximately given 
by derivatives of phase and therefore their structure corrections are unaffected 
by any constant offsets in structure phase caused by ambiguities. However, when 
fringe phase (phase delay) is the observable, ambiguities are more of a problem. 

If one attempts to define a unique value for structure phase for each value of 
(u, v) , one will find it is not always possible due to the fact that the intt)ger 
part of structure phase can depend on the path followed in the u-v plane to reach 
a given value of (u, v) , provided one requires continuity in along each (u, v) 
track. This apparent problem is really of no consequence. For purposes of cor- 
recting fringe phase in the usual VLBl applications (see Section IX and X) , it is 
not necessary to specify the integer part of structure phase since only the 
fractional part is needed to remove structure effects. After the subtraction of 
a fractional structure phase, only the corrected fringe phase has to be subjected 
to an ambiguity resolution process such as bandwidth synthesis. For this reason, 
example computations in general u-v plots will constrain structure phase to the 
range ±0.5 cycle. However, for display of results along a given u-v track, inte- 
ger cycles can be added to maintain continuity in structure phase. (If continuity 
is maintained in both corrected and uncorrected fringe phase along a u-v track, as 
is usually the case, then the procedure for structure corrections outlined above 
will implicitly enforce continuity on structure phase along the track.) 


IV. VLBI OBSERVABLES AND EFFECTIVE POSITION 

^•fith regard to geophysical/astrometric measurements, the primary observable 
in VLBl is fringe phase. As indicated in Section II, all terms but WTg, a clock 
term, and structure phase are removed from the phase (Eq. 28) to the accuracy- 
allowed by phase calibration, model restoration, etc. For simplicity, we will 
assume that all of these terms have been exactly removed so that fringe phase will 
be equal to a geometric delay term plus structure phase 
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( 36 ) 


where we have now Included the clock term t . As Indicated in Section II, t Is 
the geometric delay relative to the adopted^ref erence point ?q in the source®and 
is the structure phase given by Eq. (33). 

Three observable types are derived from fringe phase, each in a different 
way. In each case, the brightness ;transform will make a contribution, but the 
calculation of the contribution must be made for each observable using the par- 
ticular operation that is applied to fringe phase to obtain that observable. 

The three observables usually analyzed in Interferometry work are BWS delay, 
phase delay and phase-delay rate. Though it can be formulated as a finite 
difference or as a fitted phase slope, the BWS delay is approximately equal to 
the partial derivative of fringe phase with respect to frequency, (For extremely 
wide spanned bandwldths, the derivative might be an Inadequate approximation to 
the actual BWS delay. We will not analyze that problem in this report.) Phase- 
delay rate is essentially equal to the partial of fringe phase with respect to 
time, while phase delay is essentially the direct observation of fringe phase. 
Accordingly, the computation of the effect of source structure must be calculated 
as 
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for BWS delay 

(37) 

for (phase) delay rate 

(38) 

for phase delay 

(39) 


where is computed from Eq, (33). 


In general, for all three observables, the size of these contributions will 
change with baseline orientation and length so that one value for the structure 
effect cannot be quoted, even for a given baseline for a given ^ay of observations. 
Further, the contribution to BWS delay will not necessarily be the same as that 
for phase delay. Note that the structure effect depends on the reference point 
Pq, which up Co this point has not been precisely defined. A direct approach 
to removing this arbitrariness would be to assign a particular value for Pq (e.g., 
Che ordinary centroid) and then directly compute with Eqs. (37)-(39) the ^ 
structure effects for each observable as a function of hour angle, leaving Pq 
fixed. If these corrections are applied to the observables, the source would 
effectively become a point source at P^. A second. Indirect but equivalent approach 
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is more useful when detailed specification of source position is desired. In 
this approach, the effective location of the source is explicitly computed 
as a function of baseline orientation and length so that deviations from a constant 
position are directly determined. To correct the observables, a reference 
position must be chosen and the difference in effective position and reference 
position i,!.sn'nuted. The observable correction is then computed as the product of 
this differe*ice vector and the appropriate sensitivity partials relating observable 
to position. After all observations have been so corrected, the source effec- 
tively becomes a point source at the reference position. 

The effective position of a source can be determined as follows for a given 
observation. The effect of structure on our usual observables is given in Eqs. 
(37), (38) and (39). If a value of Pq is chosen that makes the structure effect 
zero for a given observable, then tbe'^observable based on Eq. (36) will come 
purely from a delay that is equal to the geometric delay produced by a point 
source at Thus the effective position will be the value of 1?^ that satisfies 

the following conditions: 


dm 


= 0 


3 <|) 


B 

at 


(}._ = 0 


m 


(41) 


(42) 


In general, the same value of Pq will not satisfy all three of these equations 
slmultaneous3.y . More exactly, one effective position can be found for the BWS 
delay and delay rate but that position will not be suitable for phase. For a 
general distribution, the next few sections derive expressions for the structure 
effects in Eqs. (37), (38) and (39) and for effective position. 


V. STRUCTURE EFFECT ON BWS DELAY 


As indicated by Eq. (37), the effect of structure on BWS delay can be obtained 
by taking the frequency partial of the structure phase in Eq. (33), which gives 



( 43 ) 
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Bosed on Eqs. (34) and (35), 


the partlals of Z become 


D cos(ioB^* l/c) 4- sin(u)Bg* l/c) j di^ (44) 

with a similar expression for Z . 

If these partlals are substituted into Eq, (43), we obtain 

At,, « • <f>„/c - B AL /c 4- B • P„/c (45) 

B s R s u s 0 

whore AL accounts for the frequency dependence of the brightness distribution and 
is given^by 
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where 
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with a similar expression for Ig. Further, the "centroid" of the resolved 
distribution is defined by the following weighted average: 




<F> , + Z/ <P> , 

Z ^ 4* Z ^ 
c s 


(48) 


for which the "centroids" of the quadrature components of the resolved distri- 
bution are given by 
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A special case is instructive. Suppose the brightness distribution is oC 
the form 


D(l^w,t) »• f(ai)DQ(P»t) (50) 


In this case, one can easily show from Eq. (46) that AL is zero so that the 
frequency dependence of the distribution has no effect on the BWS delay. (This 
result can also be obtained by considering the original brightness transform in 
Eq. (29). Note that a change in frequency for D would only change the multipli- 
cative amplitude factor f(w) and would therefore have no effect on phase.) Thus, 
if the brightness distribution retains the same shape but only changes in ampli- 
tude as a function of frequency, the BWS delay is the same as it would have been 
had the distribution had no frequency dependence. In general* since most sources 
change shape gradually as a function of frequency, the term might be very 
small for many sources. More work with actual dlatributions is required to 
assess this point quantitatively. To begin to estimate this effect, Section Xll 
develops a model for a frequency-dependent double-point source. 


VI. STRUCTlfRE EFFECT ON PHASE-DELAY RATE 


The derivation of the structure effect on delay rate observable closely 
follows that for the BWS delay. In analogy with the BWS delay case, the partial 
with respect to time in Eq. (38) leads to the result 
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Based on Eq. (35) , the partial of Z§ v/ill he given by 
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and a similar expression for 3Z^/3t. These partials can be substituted into 
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Eq. (51) to give 
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(53) 


where the centroid <p)^ is given by Eq. (48), and where ALt accounts for time 
dependence of the brightness distribution and is given by 


for which 
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(55) 


with a similar expression for H . 

For most sources, the brightness distribution will change very little 
over an observation (y 3 minutes) so that for these sources the term AT.^ can be 
neglected. Further, in analogy with the derivation for the BWS delay, one can 
consider a source for which the total power changes as a function of time but 
the shape remains the same; 


D(F,w,t) « f(t)DQ(l^u) (56) 

For such a source, one can easily show the effect of the time dependence on 
effective position is zero (ALj. « 0 in Eq. 53). Thus, although such a model is 
generally unrealistic, it does show that, if the time dependence comes from 
total power change rather than shape changes, the effect on BWS delay rate will 
be zero. 


VII. EFFECTIVE POSITION FOR BWS DELAY AND DELAY RATE 


In general multiparameter fits, BWS delay and delay rate are usually 
analyzed together, so it makes sense to specify an effective position that 
is suitable for both. As explained in Section IV, conditions on effective 
position are obtained by setting Eqs. (45) and (53) equal to zero and solving 
for i^Q, which gives 




(57) 
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fot BWS delay and 


1 • Pq * i • <i*>j^ + (58) 

for delay rate, where t is a vmlt vector In the direction of Bg. The first term 
in each of these two equations is the "centroid" of the resolved distribution, 
while the second term comes from the frequency dependence and tiii\e dependence of 
the brightness distribution, respectively. As long as Bg is not parallel to Bg, 
which is nearly always the case, the conditions in Eqs, (57) and (58) unambig- 
uously provide the effective position. However, in the general case, this 
approach vfill produce coordinate distances specified along the nonorthogonal axes 
given by B_ and B^g. To obtain coordinates along the (B, y) axes, one will have 
to apply the transformation demonstrated graphically in Figure 1. As oho’i'm, the 
effective posit^^on is given by the intersection of two lines that orthogonally 
intersect the (Bg and fg) axes at the coordinate values in Eqs. (57) and (58), 
respectively. We will not present tjie trivial linear equations that express 
(8, y) coordinates in terms of (tg, '^g) coordinates since it adds little insight 
and much algebraic litter. 


For many sources, it will probably be possible to make an approximation 
that greatly simplifies this calculation of effective position. With sources 
for which the time and frequency variability are negligible (ALt. » 
effective position will be given by 


Pq - (59) 

This can be easily seen by considering Eqs. (57) and (58) with Ab^ and AL^^ set 
equal to zero. Thus, in this important case, the effective position becomes the 
centroid of the resolved distribution. 


A special case is worth giving as an example. For sources that are compact 
relative to fringe spacing (i.e., Bg/1 0), one can easily show that the 
centroid of the resolved distribution becomes 
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Further, one can show AL » AL * 0 when B A *>■ 0, Thus, if the source is very 
compact, the effective position becomes thi ordinary centroid of the distribution. 


Vm. EFFECTIVE POSITION FOR PHASE DELAY 


The effective position for the phase observable can be obtained from Eqs» 
(S3) and (42), which gives the condition 



where Z and Z are given by Eqs, (34) and (35). Once the brightness dlstrlbu- 
tion^is'^aDeclfled, Z and Z can be readily determined, thereby yielding a value 
Cor B • K. Thus, tor the brightness transform to make no contribution to phase 
delayf the^component of effective position along must be made egual to the 
expression in Eq. (61). In general, the effective position along B will not 
be the same for the phase-delay observable as it is for the BWS delay observable. 
As in the case of BWS delay, however, the phase delay places no requirement on 
the position component in the direction orthogonal to B^. The next section will 
discuss the correction of phase delay for structure effects. 


As was done in the last section for Che effective position for BWS delay 
and delay rate, one can consider the special case of zero baseline. For zero 
baseline one can easily show that Eq. (61) becomes 


A 



■^0 • \ • <’^>0 


(62) 


where is given by Eq. (60). 'flius, as for the other observables, the 
ordinary centroid will again give the effective position when source extent is 
much smaller than the fringe spacing. 


IX. STRUCTURE CORRECTIONS AND THE ORIGIN OF STRUCaURE COORDINATES 


Mien source structure measurements are made, the origin of the r.a. and 
declination coordinates (structure coordinates) over which Che brightness is 
distributed is not accurately known relative to absolute coordinates, even though 
the relative location of features in the distribution can be accurately specified. 
This section shows that, if the brightness distribution is "perfectly" determined 
in this relative sense, the absolute location of the origin of structure 
coordinates can be determined, in principle, to the accuracy allowed by the 
errors in measured delay and delay rate, as propagated through a standard multi- 
parameter fit to determine baselines, source locations, etc. An important 


limitation will be the ercor in aligning the S-band and X~band distributions 
and the resulting delay error due to S/X calibration of charged particle delays. 
In the process of explaining the determination of the origin of structure 
coordinates, this section also shows how brightness distribution measurements 
can be used to correct for the effect of structure on delay and delay rate. We 
will assume that the extent of the source is sufficiently limited and that the 
origin of structure coordinates is fairly close to the center of the source so 
that second order effects for partials w.r.t. position, etc., can be neglected. 
Further, ic will be assumed that BWS delay and phase delay rate are the observ-^ 
ables. The case of phase delay is discussed at the end of the section. 


Since origin determination is closely associated with position determination, 
we will choose to use the concept of effective position to justify the procedure 
for absolute determination of the origin of structure coordinates. Although 
we will do so in this report, it will be shown that effective position does not 
necessarily have to be determined explicitly in practice. 


We will first consider the case in which charged particle delays are 
negligible. As shown in Section IV, the fringes can be formulated so that the 
brightness transform makes no contribution to the observables of BWS delay and 
delay rate. This is accomplished by setting the reference position equal to 
the effective position discussed in Section VII. When this is done, the resulting 
observables of BWS delay end delay rate, which are obtained through the differ- 
entiation of Eq. (36) with respect to frequency and time, will be given theoreti- 
cally by 


■'bus “ ''s ■'c 


• • • 

(}i g C 


(63) 


(64) 


where tg is the standard geometric delay for the effective position . Explicit 
structure contributions (8(j)jj/3ti) anc^ D(j)B/3t) have disappeared as a consequence 
of setting the reference position Fq equal to the effective position. The most 
Important aspect of this result is that both observables are the observables 
that would be generated by a point source at the effective position. 


At this point, the absolute location of the effective position is unknown 
but, given the brightness distribution, the location of effective position can be 
calculated accurately relative to the origin of structure coordinates (or any 
other point in structure coordinates), as described in Section VII. If one wishes 
to apply a differential correction to the observables that will shift the effective 
position to another point in the source, he can do so by computing and applying 
a correction equal to the product of the vector difference in position (computed 
in structure coordinates) and the sensitivity partials relating observable to 
position. When such differential corrections are applied to the total observables 
in Eqs. (63) and (64), the new geometric delay and geometric delay rate will 
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pertain to a point source ploced at the new reference position. The absolute 
coordinates of the new point will still be unknown while its position relative 
to structure coordinates will be accurately known. 


The absolute coordinates of the new reference position can be detemtlned 
as follows. Typically, a measurement session will Involve several observations 
of a given source, for which the effective position changes from observation to 
observation, For each observation, the observables will be corrected for the 
displacement (the solution to Eqs, 57 and 58) of effective position from the new 
reference position. We will assume for each source that the new reference posi- 
tion Is placed at the structure-coordinate origin, although another point could 
be chosen. Wlien the observables are passed through the usual multiparameter 
fit for source location, baseline, etc, , the solve-for source position for each 
corrected source will be the absolute location of its origin of structurs 
coordinates. The brightness distribution for each source has how been placed 
absolutely on the celestial sphere to the accuracy allowed by the errors In the 
observables as propagated through the multiparameter fit and by the errors in 
the measured brightness distribution. If the data covers a long time span so 
that the brightness distribution of a given source chlinges significantly and must 
be remeasured, then one must face the problem of determining the origins of two 
or more sets of structure coordinates for one source. Two choices present them- 
selves. One can either solve for a separate position (origin) each time a new 
distribution Is Introduced, or all origins for a source can be forced into 
agreement by prominent feature alignment. The choice of an approach depends on 
experimental goals and the nature of the source. 

The preceding discussion was presented to demonstrate through position 
analysis how the origin of structure coordinates can be absolutely determined. 

As one might guess, it is not necessary to compute effective position explicitly 
but only structure delay (and delay rate) relative to the origin of structure 
coordinates (Eqs. 43 and 51 with - 0). Following the procedure outlined 
above, one would need the sensitivity relation given by 


At ^ -Bg • ht/c (65) 


where AS is the shift In position for the source. Note that If the effective 
position described in Section VII is substituted in_^place of Aa as prescribed, 
the resulting delay will be given by Eq. (45) with Fq ■ 0, by design of the 
effective position. A similar argument can be made ror delay rate. Thus, if 
one computes the structure delay (and delay rate) directly from the derivatives 
of structure phase (Eq. 33) with * 0, the resulting delay (delay rate) is the 
correction needed to make the observable equal to that generated by a point 
source at the origin of structure coordinates. 


When charged particle delays are considered, the procedure becomes somewhat 
complicated, even thoxigh the basic concept^ are the same. If observations were 
made at one frequency and outside calibrations could supply accurate values for 


19 



charged particle delays, the procedure would be the same as that explained 
above. However, if we assume dual frequency (S/X) calibrations of those 
delays, we mjst complicate the procedure with two brightness distributions, 
one from each frequency, Tlie main problem is that each distribution will have 
its own unknown origin unless steps have somehow been taken to align them or to 
determine their relative placement. At this time, a method has not been veri- 
fied that will accurately align origins in the general case. We cannot use the 
preceding approach separately for S-band and X-band since the charged particle 
delays will corrupt the determination of absolute position, Though it may not 
be very satisfactory, one solution is to use the two measured brightness dis- 
tributions to align prominent features via a translation of the S-band distri- 
bution, The accuracy requirements on such a translation will be very stringent 
but, due to the nature of S/X calibrations, it will be reduced by a factor of 
13 relative to the standard error requirement. For example, a 1-cm error in delay 
on a 6000-kin baseline would be produced by a misalignment error of l3*0.0003*Vcm * 
1 cm >• 0V004, Determination of whether enough sources can be treated by the 
"prominent-feature" approach will require a study of actual S/X structures, 

'f'nure are apparently a number of sources for which this would be a feasible 
approach (e,g,, Ref, 2), 


Once the origin of the distribution at S-band has been made consistent with 
that at X-band, whatever the approach, the next step would be to correct each 
obsarvable for the displacement of effective position from the origin of struc- 
ture coordinates, as described above, The S-band and X-band observables would, 
of course, be corrected according to their respective brightness distributions. 

Due to differences in structure, the two bands will typically have different 
corrections. This procedure would, in effect, make the source a point source 
at both S-band and X-band, with a location equal to the (presently unknown) 
position of the origin of structure coordinates at X-band. After the completion 
of this procedure for all sources, the next step would be to perform S/X cali- 
bration for each observable type by means of the standard linear combination of 
S- and X-band observables. Such a calibration would now be valid since the S- 
and X-band positions have now been made coincident for all observations, pro- 
vided the aforementioned alignment of origins has been accurately carried out. 

The corrected observables (the output of he S/X linear combinations) would then 
be simultaneously subjected to a multiparameter fit for baseline, source position, 
etc. The resulting solution gives for each source the absolute location of the 
(S/X common) origin of structure coordinates. As a result of this procedure, 
the S- and X-band distributions have both been accurately placed in an absolute 
sense on the celestial sphere for all sources. 


The preceding discussion assumed that BWS delay and delay rate were the 
observables. If absolute phase delay is the observable, the concepts are 
largely the same but the equations change. In analogy with the other observables, 
the phase delay will be corrected by (|>g in Eq, (33) with Fq 0. Unlike the 
structure corrections for the other observables, structure phase is afflicted 
with integer-cycle ambiguities. The method for handling these ambiguities 
consists of the following steps. First, the fractional part of Is computed 
through Eq. (33) and subtracted from fringe phase, which itself still possesses 
an absolute ambiguity. This ambiguity in corrected fringe phase would 
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then be removed by the BWS technique in which an iterative process successively 
removes the ambiguities in (structure-corrected) BWS delays for ever-wlder 
channel separations until RF phase la reached. For the BWS procedure to succeed, 
the structure-corrected BWS delays and phase delay must all be equal to the 
same delay, except for ambiguities and measurement uncertainties. Therefore, 
to equalize the geometric component, structure effects in BWS delays must be 
computed and removed relative to the same assigned source position as is phase, 
(The computation of structure effects in Section IV approximates a finite 
difference with a derivative. For very widely spaced channels, that approxima- 
tion might fall. It is beyond the scope of this report to treat the differential 
analysis.) After absolute removal of integer-cycle ambiguities through the BWS 
technique, the resulting phase delay will be equal to the geometric delay that 
would be generated by a point source at the origin of structure coordinates. 

S/X corrections also closely follow the preceding discussion for BWS delay and 
delay rate. 


Associated with the preceding consideration of structure effects is the 
issue of assigning a nominal location to a source. Once milllarcsecond accuracies 
are reached, two numbers will no longer be sufficient to specify the location 
of a complicated source. That is, the effective location will vary as a function 
of baseline length and orientation. As suggested above, source geometry has been 
completely specified once one has been supplied the brightness distribution and 
the absolute location of the origin of structure coordinates. Since the origin 
can be arbitrarily shifted to any place in the brightness distribution, one 
should attempt to make the two required numbers specifying the origin have as 
much utility as possible. One can make the point that the "compact component" 
should be made the origin since long-baseline solutions will apply to the compact 
component and will give the most accurate positions that are to be compared. 
However, in the most accurate measurements, the "compact component” position can 
be strongly dependent on the geometry of resolution and therefore only has 
meaning In terms of a specific geometry. One attractive possibility is the 
ordinary centroid, since it gives, the Average location in the conventional sense 
and is the correct effective position for the zero-baseline case. Further, it 
requires no specification of baseline orientation to give it meaning. It is 
definitely better than an arbitrary assignment, since effective positions computed 
in structure coordinates would then automatically be given in terms of deviations 
from the zero-baseline case. 


X. DIFFERENTIAL POSITION MEASUREMENTS 


One application of differential VLBI is the measurement of the relative 
positions of two natural sources. By forming the difference of the observables 
for closely spaced sources, one can remove many errors, thereby allowing very 
accurate determination of differential positions. The resulting differential 
observables are fit with a delay model possessing only differential right ascen 
sion and declination as solve-for parameters and perhaps a constant "clock" 
term (for the phase-delay case) . 
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As Indicated earlier, source structure effects can be quite large and must 
be corrected in some cases. The correction procedure would be like that described 
in the last section, but now two distributions would require corrections, one for 
each source. Thus, two structure corrections would be applied to each differ- 
ential observable, where aach correction is computed relative to some assigned 
reference position. In effect, this turns each source into a point source located 
at its respective reference position. When all observations have been corrected 
and then subjected to a multiparameter fit, the resulting differential coordinates 
will refer to the difference in the two adopted reference positions. 


At this point we are no longer dependent on a particular observation and 
geometry for a specification of the relative location of the sources. Since we 
now have a complete description (i.e., the two distributions along with an ac- 
curate measure of the separation of their origins), we can now select whatever 
measure of separation we choose. This report will not propose a "best" definition 
of that measure since more consideration of the problem is required. One possi- 
bility would be to use the "compact-component centroid" for a reference geometry 
of resolution, where the same geometry would be applied to successive experiments. 
In this way, changes in resolution geometry would be eliminated as a source of 
error. However, changes in the compact components could still contribute sub- 
stantial errors to the measured change in differential effective position between 
experiments. Such effects would have to be corrected through comparisons of the 
full distributions, perhaps by prominent compact-feature alignment between the 
brightness distributions of successive experiments. 


A few words are necessary to explain the treatment of integer-cycle ambigui- 
ties in structure phase and fringe phase in the case of differential measurements. 
To obtain the final observable, it is only necessary to compute the fractional 
part of (f)j^ in Eq. (33) and subtract it from fringe phase. (However, one can 
restore integer- cycles to structure phase in a relative sense by requiring con- 
tinuity along a u-v track.) The ambiguities still present in corrected fringe 
phase can be removed in one of three ways, depending on the type of experiment. 

In single-channel measurements with poor a priori values for baseline, etc., 
ambiguities can only be removed in a relative sense by requiring continuity be- 
tween observations. The unknown overall bias created by unresolved ambiguities 
would be modeled and solved-for in parameter estimation. In slngle-^channel ex- 
periments with sufficiently accurate a priori information, ambiguities can be 
removed abso.!utely by comparing each measured value with differential phase 
predicted on the basis of a priori information. In multiple-channel experiments 
with suitably placed channels and phase calibration, ambiguities can be removed 
through the BWS procedure. For the last two absolute schemes, there would be no 
need to solve for a bias term provided adequate differential cancellation of 
Instrumental effects, etc., has occurred. 


Due to small observable errors and geometric considerations, differential 
measurements will provide a good direct test of the theory of structure effects. 
Measurements of delay, for example, are^directly proportional to the displacement 
of effective position in the direction for each observations. Because of 
this, delay residuals and their reduction through the application of structure 
corrections can readily be Interpreted geometrically in terms of effective position. 
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XI. THE SYMMETRIC BRIGHTNESS DISTRIBUTION 


For a symmetric distribution, one can easily show that the structure phase 
is some value of mr, provided the center of symmetry is taken as the origin of 
structure coordinates. By definition the distribution is symmetric if 


D(S,Y»w,t) » D(-S,“Y,w,t) (66) 

If in Eq. (24) we place Yq the origin of structure coordinates 
(Sq “ Yq “ 0) and take the complex conjugate, we obtain 

R*(u,v,ai,t) D(S,Y,w,t)e 2iri(uS+vY) 

By a change of integration variables given by 

S' » -3 (68) 

y' ® -y (69) 


and by the use of the symmetry relation in Eq. (66), we obtain 


R*(u,v,..,0 dSMy' 


(70) 


which is the complex conjugate of Eq. (67). Thus, when the reference position is 
placed at the center of symmetry of a symmetric source, the brightness transform 
R is equal to its complex conjugate. This shows that, for a symmetric distribu- 
tion, the transform is real and the structure phase is therefore equal to mr. 

In the limit of zero baseline length (u = v = 0) , the brightness transform becomes 
the total power of the source, which is a positive real number. In this limit, 
the transform phase will be defined to be zero. Since the transform phase of 
a symmetric distribution can only be equal to the aforementioned discrete 
values of mr, it is clear that if phase changes, it must change in jumps. However, 
for a well-behaved brightness distribution, the brightness transform will be a 
continuous function along a (u-v) track and consequently the transform phase can 
only change at those u-v values for which the transform (fringe) amplitude is 
zero. The transform phase will be constant along tracks away from those u-v values 
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and, If it changes, the change must occur in the form of a discontinuous jump 
at a null amplitude point. We will see an e;;ample of such behavior in the follow- 
ing section. Since BWS delay and delay rate are derivatives of phase, the effect 
of a symmetric distribution on these observables will be zero except at points 
of zero amplitude, where the effect can be Infinite (6-function response). For 
the more realistic case of a nearly symmetric distribution, the structure effect 
on BWS delay can be quite large near (u,v) points of nearly null amplitudes and 
quite small away from such points. It follows from the preceding discussion 
that, if the brightness transform of a symmetric distribution possesses no 
finite (u,v) points of null amplitude, then the transform phase is zero for all 
(u,v) . 


XII. THE DOUBLE-POINT SOURCE 


The next example will be a time-invariant source consisting of two close 
point sources whose strengths and spectral indices can be varied, The brightness 
distribution will be described by 


D(^,w) 



6 (f - f + 



6(P + 


(71) 


where we have chosen the midpoint between sources as the origin of structure 
coordinates , and where k„ is the spectral index of source A. The factors 
can be used to vary the ^ strength of the sources, while Wq is a normalizing 
reference frequency. For this distribution the brightness^transform in Eq. (29) 
becomes 
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Fringe visibility can be easily computed from this expression by taking the 
RSS of the real and Imaginary parts and dividing by the total flux, g^^ + g2> which 
gives 


V =Cv^+(l-v^) cos^ (wB • 
s m m s 



(75) 
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where the minimum fringe visibility is given by 


V 

m 


" ^2 
8^ + 82 


( 76 ) 


(Note that v can be either positive or negative so that, strictly speaking, the 
absolute value of v is the minimum fringe visibility.) Further, one can easily 
solve for structure*" phase to obtain 


(jijj - -tan“^ jvj^^ tan (wB^ • P^^/c)] + wB^ • Pq/c 


>] 


(77) 


It is convenient to rewrite this equation in the form 
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*• “tan |v tan 
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^v^^ tan (0/2) j + 6/2 


(78) 


where 


0 » 2irSg/A 


(79) 
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The quantity Sg is the component of the vector 6S “ 2P^ containing the two 
sources projected onto the baseline vector and is given by 


S 


B 



( 80 ) 


The quantity A is the Interferometer resolution (spacing of fringes on the sky) 
and is given by 


A 


- 2it c 
~ W B 

s 


(81) 


(82) 


We have chosen this representation because it conveniently makes 0 (in cycles) 
equal to the projected separation of the source divided by the interferometer 
resolution. Further, in Eq. (78) we have set so that the phase that will 

be plotted will be referenced to the location of point 1. With this choice of 
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riSferonce point, the phase oscillates about zero, remaining between +7 t/ 2, pro- 
vided source 1 is the stronger source. Figure 2 shows the behavior of <|)g for 
various values of the strength ratio 82 / 81 * This phase behavior can be easily 
understood by plotting the phasors in Eq, (73) in the complex plane with Fq " F,, 
as schematically shown in Figure 3. For 82/61 ^*^ (1*^*» equal strengths)t 

the phase forms a sawtooth pattern as a function of S^/A with jumps of ir whenever 
the projected source separation is an integer multiple of one-half of the fringe 
spacing. In terms of phasors in Fig. 3a, such jumps occur whenever the second 
phasor has a phase of nrr so that the sum amplitude is nearly zero and the sum 
vector moves rapidly from the first to the fourth quadrant as a function of the 
phase of the source - 2 phasor. When the strength of source 2 is very weak 
compared to source 1 ( 82 / 6 ]^ < 1 )» the phase oscillates by a small amount about 

zero, as shown for 82/61 ®* 0*1 in Fig. 2 and illustrated with phasors in Fig. 3b. 


Effective position for phase-delay, relative to the origin of structure 
coordinates, can easily be obtained by setting Eq. (77) equal to zero and solving 
for Pq: 


where S is the separation of^the two sources (S « 2|P. j). Thus the component 
of effective position along B is a fairly complicated function of resolution 
geometry and source characteristics. Numerical examples of effective position 
are given in Fig. 4 under the assumptions that the baseline vector D points along 
the vector connecting the sources (Sg**S) and that structure phase is^a continuous 
function of S/A that equals zero when the baseline length equals zero. As in- 
dicated in Eq. (83), the results are given in Fig. 4 in units of source separation 
so tliat deviations from the position of the stronger source (//I) are to be com- 
pared source separation sho\m in the figure. Note that the effective 

position 'Uiu’les from source 1 by less than one-half of the source separation 
and decreases as the strength of source 2 decreases. The decrease in the oscil- 
lation of effective position with increasing S/A is a result of increasing 
sensitivity of phase to position change (e.g., 8«})/9d becomes larger as baseline 
length increases). Though the effective positions shown here for phase delay are 
not as constant as one would like, they are much more stable than the effective 
positions that will be derived for BWS delay later in this section. 

The effect of structure on BWS delay can be obtained by differentiating Eq. 
(77), with respect to to, which gives 



9to 


+ T 


+ B 




(84) 


where 
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Figure 2. Structure Phase* for a Double-Point Source 
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PHASE DELAY AT 8.5 GHz, psec 
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Figure 4. Effective Position for a Double-Point Source: Phase Observable 
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— — ^ tan (2irBg • Pj^/X) 


■^D = 2u 


(85) 


1 -f h‘ 


t:_ s !L_ ^ h sec (2 itB^ • fJX) 

4 C Si. 


p 1 + h' 


( 86 ) 


»> -V 


h 5 tan (2rrBg • V^/\) 


(87) 


The first term comes from the frequency dependence of the brightness distri- 
bution while the*^second term tp comes from the frequency dependence of u, v 
in the transform or, in other words, from the frequency dependence of the sky 
fringes. Note that, as one would expect, the first term is zero when « k 2 « 
(t'Jhen kj^ “ k 2 » the shape of the distribution does not change with frequency so 
that the BWS aelay effect due to frequency dependence is zero, as explained in 
Section V.) Further, the second term Tp is zero if « g 2 > except at singular 
points. (When g.|^ » g 2 » the distribution is symmetric about the origin of struc- 
ture coordinates’. As^explained in Section XI, such symmetry leads to zero effect 
on BWS delay, except for possible singular points.) 


The effective position associated with the BWS delay can be obtained by 
finding the value of ? 0 ->.that makes the delay in Eq. (84) zero. The value of tlie 
component of along that will do this given by 


1-fh 


X g-v^) 

B 4r 


tan (2 ttB 


P-/1) + v^ 

X ms 


^ 2 ^ 
P, sec (2 ttB 
1 ' s 


Pj^/A) 


( 88 ) 


More will be said about this result for effective position after the delay rate 
has been analyzed below. 

Before proceeding to the delay rate observable, numerical examples will be 
presented for the two delay terms in Eq. (84) that give the effect of double-point 
structure on BWS delay. The first term in Eq. (84), the one due to the frequency 


r— 4c»=e-sc!sr.'wrt' 
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dependence of the point strengths, can be 'rewritten as 


■'d " 


(k^ - k,) (1 - 


2(13 


1 + v^“ ton" 
tn 


( 0 / 2 ) 


( 89 ) 


where 0 was defined above as the ratio of projected source separation and 
resolution. Since this funcfcion is periodic and antisymmetric about 0 “ 0, one 
plot of Tjj for 0 < e/2 < it/ 2 suffices for all 0, Figure 5 plots Tj. as a function 
of 0 for various values of relative source strength gn/s^* with the spectral 
index difference equal (k, - k^^) to 1.0. Since the effect is proportional to 
*^2 ~ ^1’ delay for other values of ^2 - can be obtained by multiplying the 
results from Fig. 5 by k2 ~ k^^. The delay effect is zero for 0 ■ 0 (or 2mr) and 
passes through a maximum toward r (or (2n + Dr). Note that the largest values 
(e.g. , > 0.2 nsec) are obtained when the source strengths are nearly equal (0.95 
< g^/Si I'O) and the spectral index difference is large (k„ - kj^ > 1.0). An 
important quantity is the duration of the delay "pulse" as a^function of time 
(through 0) on an intercontinental baseline. For a 6000-kn» baseline observing a 
source with separation of OVOOl at X-band, 0 can vary by no more than 0.005 
cycles in 1 minute. Thus, in this case, only a fraction of the "pulse" width 
will be traversed in a 3-minute observation, when 82/8^^ 0.95. 


A general expression can be easily computed for the maximum value of Tw as 
a function of 0. Such a maximum occurs whenever 


tan(0/2) " i 

m 


at which point the maximum of Tjj Is 


D 


max 



(1 - V ^) 
2 m 

V 

m 


(90) 


(91) 


This function is plotted in Fig. 6 at 8.5 GHz as a function of the spectral index 
difference (k^ - k2) and of the strength ratio . Note that the delay 

effect can become infinite as the strengths converge (g2 gj)* The reason for 
this singularity can be seen in Fig. 3a. For g2 = gx and 6 z (2n + 1)7T, the 
brightness transform is the sum of two phasors of nearly equal length pointing 
in nearly opposite directions. Consequently, the amplitude of the sum phasor is 
very small compared to the amplitude of the component phasors, and a slight 
change in the amplitude of the component phasors (i.e., of the point strengths 
gj^ and go) as a function of frequency (as in the partial 3^g/3co) can cause large 
changes In the phase of the sum phasor. As we shall see, similar singular be- 
havior is exhibited by the second term t in Eq. (84) as a consequence of the 
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Maximum BWS Delay* Due to Frequency Dependence of a Doutle-Foint Source 


phase chongea In the component phasors. The special case suggests chat if a 
given source exhibits very small fringe amplitude relative to the amplitude for 
other (u,v) points, then that source should be investigated for undesirable 
behavior in structure phase at that amplitude minimum. For the particular case 
of a double point source with ^2 - kx » 1.0, the maximum value of the structure 
delay term tjj reaches 0.18 nsec (5.5 cm) at x-band when g2/si " 0.95 (or 
visibility •• 0.026), 


To sununarlze the results for this delay term, the frequency dependence of 
structure can theoretically cause significant delay errors in the case of special 
double-point sources, if they arc not corrected. However, these worst-case 
sources would probably be unusual since they would have to possess both nearly 
equal point strengths and a large spectral index difference. Given adequate 
structure measurements as a function of frequency, corrections can be applied to 
the delay observables to remove this effect adequately. 


For the same example of a double point source, the second structure effect 
(t ) on BWS delay is caused by the frequency dependence of the sky fringes, given 
by^Eq. (86). This term can be rewritten as 


T 


I 

P 


ZhL 8/2 ^ ^ 

cos^(0/2) + v^^sin^(0/2) ^ ^ 


(92) 


where v and 0 are defined above and where the second term has been included 
in order to reference the delay to the brightness centroid given by 



(93) 


Though such a subtraction is not necessary, it removes a distracting linear 
drift (with respect to 6) from the delay and comes close to referencing the 
structure delay to an "average" or "typical" source position. This is equivalent 
to assuming that the last term in Eq, (84) has been combined with t where we 
have made the definition Pq *“ P^. Figure 7 plots t' over a range ^ of 0 at 
8.5 GHz for few values of '^the strength ratio g2/g^? Relative maxima occur near 
0 » (2n4*l)7r, with the maximum values increasing with Increasing 0. A fairly good 
approximation for these maxima can be obtained by setting 0 « (2n + 1)tt, which 
gives 


X 


I 


tt(1 - V ) 
m 

wv 
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n 

2u 


max 


(94) 
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where 


W" (t) " n - 0, + 1, + 2, ... (95) 

n 

Thus the maxima increase in magnitude as the source strengths converge (for 
the reasons discussed above) and as the separation/resolution ratio increases. 


Figure 8 plots t' 


max 


at 8.5 GHz as a function of S^/A and 


- 0.95 and separation/resolution « 0.5, the maximum delay effect is about 
I nsec. For a wide range of strength ratio and S^/A, the effect is greater than 
30 psec (1 cm). Thus, it appears that, for many sources that are nearly double 
point, this structure effect can be very significant and must be corrected if 
accuracies at the centimeter level are required. The problem at ids worst levels 
would be alleviated, of course, by avoiding observations in the regions of rela- 
tively small fringe amplitude. 


Given the preceding analysis for a double-point source, one can compare the 
magnitude of the delay effect due to the frequency dependence of the distribution 
(Xjj) with that due to the frequency dependence of the sky fringes (t ) . A 
comparison of Eqs. (91) and (94) indicates that t is about an order^of magnitude 
larger than x^. This result suggests that one meisurement of the brightness 
distribution at the observing frequency will be adequate to remove most of the 
structure delay and that a second brightness measurement to determine frequency 
dependence will be necessary only for the most accurate applications, if at all. 


In the same example of a double-point source, the effect of structure on 
delay rate can be obtained by taking the partial of (ji^in Eq. (77) with respect 
to time: 



V 

m 


1 + h^ 



sec ^(2Trl^ . + Bg • Pq/c 


(96) 


where we have assumed the brightness distribution is independent of time. Since 
the delay rate observable is much less important than the BWS delay as a conse- 
quence of larger measurement errors, we will not present numerical examples of 
delay rate. However, comparison of Eqs. (86) and^(96) will reveal that, with 
respect to occurrence of maxima, the behavior of x, will be similar to the BWS 
delay. \^hen converted to equivalent baseline error, the sizes of these maxima 
are of the same order of magnitude as the maxima in x in Fig. 7. 
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FLUX 1 


Figure 8. Maximum BWS Delay* for Frequency-Independent 
Double-Point Source 
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By setting the derivative in Eq. (96) equ^l to zero and solving for Pq, 
we obtain another condition on a component of P^; 
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P, 
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sec^ ( 27 rB 


1+h' 


—i B 


s 


( 97 ) 


Both BWS delay and delay rate are normally used together, in which case Eq. (88) 
and Eq. (p7) can be used to determine uniquely the effective position P^. As 
explained In preceding sections for a general distribution, that procedure is 
usually somewhat involved if the frequency dependence of the brightness distri- 
bution is not negligible (Xj^ * 0). l^en such dependence is negligible, the 
effective position has been shown to be equal to the centroid of the resolved 
distribution. In the present example, this centroid is easily obtained by 
inspecting Eq. (97) and the second term in Eq. (88) which, along with Eq. (87), 
gives 



V f, 

m 1 

cos^(27r§ • /A) + V ^ sin^ (2irf • f~/X) 


(98) 


As emphasized after Eqt (71), this position calculation is referenced to the 
midpoint between point sources. As one would expect, the centroid lies on the 
line connecting the two sources, but is a rather complicated function of the 
source strengths and resolution geometry. Three special cases are of interest. 
For various values of the minimum fringe visibility, it is easily shown that 


<?> 

<^> 

<f> 


0 where g^ » g 2 [v^ = 0 and 0 (2m + 1 )it] (99) 

when g^ » g 2 (v^ = 1) (100) 

when g^ « g 2 (v^ = -1) (101) 

^ when S = 0 (102) 

c s 


as one would guess. 


To compute numerical examples of the effective position, it is convenient 
to rewrite the components of Eq. (98) in the form 



1 


(103a) 


<g> . % 

® ^ cos^(e/2) + V^^sln^ (0/2) 

<Y> “ 0 (103b) 


where S and 0 have been defined above. Without loss of generality, we have 
placed along the 3-axis and have expressed results in units of S, the distance 
between sources. In this form, the result depends only on the strength ratio 
( 62 / 61 ) ^t^d on the ratio of projected separation to resolution (0). Figure 9 
plots'*'ef fective position as a function of 0 for a few values of the strength 
ratio. As can be seen, the effective position is a periodic function of separa- 
tion/resolution that goes through a maximum whenever separation/resolution 
« (2n + l)/2 cycles. The function becomes singular at the same points if the 
point strengths are equal, for the reasons discussed earlier. Even when 82 / 6 ]^ 

« 0.5, the effective position can move outside the region between sources 
to the extent that it is separated from the stronger (closer) source by a distance 
equal to the source separation. Thus, the movement of the effective position is 
generally much more damaging in the case of BWS delay and delay rate than it is 
in the case of phase-delay (Fig. 4), This result suggests that, with respect to 
structure effect, the phase delay will tend to be superior to the BWS delay. 


Due to the importance of null points and minimum amplitudes, it is useful 
to demonstrate their location in the u-v plane for a double-point source. As 
can be readily discerned from Eq. (75), minimum amplitude will occur on lines 
for which ^„/X satisfies 

P 


fs ; 2 n + 1 
X 1 ■ 2S 


(104) 


where S is the spacing in radians between point sources (S = 2|Pj^[). Figure 10 
schematically plots these lines of minimum amplitude in the u-v plane. As in- 
dicated by Eq. (104), the lines are perpendicular to and have a spacing of 
S”^, with the first lines separated from the origin by (2S)“^. It is clear 
from this result that, for a double point source, a long u-v track will usually 
cross a line of minlnmm amplitude if /X is greater than (2S)“^. Thus for nearly 
equal point sources, it will not be unusual in such cases to apply fairly large 
structure corrections to BWS delay for some observations. 
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Figure 9. Effective Position for a Frequency-Independent Double-Point 
Source: BWS Delay and Phase-Delay Rate 
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Figure 10. Lines of Minimum Fringe Amplitude for 



XIII, THE N-POINT SOURCE 


This section analyzes the structural effects of a frequency-independent 
source consisting of N-point sources of arbitrary strength and position. The 
special case of a triple point source is consUered in some detail. 


The brightness distribution for an N-point source can be written as 


N 

D(t) - 2 Bj, - V 

k*l 

where gj^ and are the strength and position of point source k. For this 
distribution, the brightness transform in Eq. (29) is given by 


K.B • ¥ /c * -i,.B . P./c 

„ so X s k 

R - e 

k=l 


(106) 


for which the structure phase becomes 


<|.jj = tan‘^(-Zg/Z^) + o)Bg • Pq/c 


(107) 


where 




(108) 




(109) 


®k = 


( 110 ) 
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The fringe visibility becomes 



(111) 


The BWS delay due to structure is obtained through the use of Eq. (45) with 
» 0, which gives 


m “ 

ATg « - • tjc (112) 


where the effective position is given by 


c s k 

Thus the effective position is equal to a linear combination of the point-source 
position vectors where the coefficients are complicated functions of source 
and resolution geometry. For the special case of N colinear xioint sources, this 
result shows that the effective position lies on the line joining the sources, as 
one would expect. The special case of two point sources can be used as a check 
by showing Eq. (113) reduces to Eq. (98) for N = 2. 


Unlike sources for which N is greater than 3, a triple point source can be 
characterized by a sufficiently small number of variables to allow a reasonably 
general display of results. As Indicated in Eq. (106), the brightness transform 
for a triple point source will consist of the following sum of three phasors: 


R * 1 + g2 e”^®2 4- (114) 

where, for simplicity, we have used_^the ^trength and position of point source 1 
as reference by setting gj^ « 1 and Pq » Without loss of generality, we will 

assume point source 1 is the strongest source. Thus, four variables can be used 
to parametrize a triple point source: g 2 , g^* 02> ®3* Us schematically 

illustrates the formation of the brightness*^ transform from the three phasors. 

As the baseline changes for a given source, the directions of the phasors (i.e. , 
02 and 6_) will change while, of course, their lengths (i.e., point strengths) 
remain constant. The resulting length (amplitude) and orientation of the sum 
phasor will change as a function of 02 and 0 . This parameterization in terms 


43 


(a) 


THE STRONGEST POINT 
SOURCE IS REFERENCE 



REAL 


Figure 11, Schematic Example of (a) the Formation of the Brightness 
Transform for a Triple-Point Source and (b) Zero Fringe 
Amplitude for a Triple-Point Source 
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of (02» 63) will b6 general in that che resulcing pliaae and delay equations 
will be valid for any values of source positions and baseline. To use the 
result, one merely makes the conversion 


( 115 ) 


( 116 ) 


where Bs/1 and (Pg, P3) are the given values for the baseline vector (i.e., u, 
v) and the source positions, respectively, 


- 2" f • ‘’2 

-y 

“3 ■■ • '^3 


Since large delay effects occur where the amplitude of the brightness trans- 
form approaches zero, we will first determine where points of zero amplitude 
occur. Zero amplitude occurs whenever the phasors sum to zero as shown in Fig. 
libs For given values of point strength (g2, gg), there are unique values of 
(021 0 ^) which such closure will occur. One can easily show that those values 
are given by 


2 m 7 f m “ 0, + 1 , +2, ... ( 117 ) 


2nif n 0, + 1, +2, . . , ( 118 ) 


where the leading plus (minus) sign in the second equation is used if the leading 
minus (plus) sign is used in Eq. ( 117 ). Such null points can occur only if the 
two weakest point sources have a combined strength at least as strong as the 
strongest source (i.e., 82 + 83 i !)• Figure 12 schematically demonstrates how 
the locations of these null points can be calculated in the u-v plane. In the 
plane of the sky at the source, draw two new axes in the directions of point 
sources ( 2 , 3 ) starting at point source 1 (I.e., along P2, Fo). Then for point 
source 2, draw perpendiculars to the ¥2 axis that intersect chat axis at distances 
from the u, v origin equal to 





(m) 

, m = 0, + 1, + 2, ... 


( 119 ) 
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FLUX 2 + aux 3 > aux 1 (LARGEST) 



Figure 12. Schematic Example of the Locations of u-v Points of Zero Fringe 
Amplitude for a Triple-Point Source 




Next, draw similar lines perpendicular to the P« axis, based on ? 3 , etc. The 
intersections of these two sets of perpendiculars are the location of the null 
points in the u-v plane. As can be seen, the null points occur at discrete 
points on a regular lattice. For simplicity, the figure shows only the null 
points associated with the leading plus sign in Eq* (117), The other set of null 
points that comes from the leading minus sign is also latticelike but is generally 
not coincident with the first set. 


Tlie null -point behavior of a triple-point source is to be contrasted with 
the result for a double-point source for which the zero amplitudes fall on lines 
in the u-v plane (see Fig. 10), For a double-point source, the u-v points of 
minimum amplitude depend only on the source separation vector and not on point 
strengths, On the other hand, for a triple-point source (and for N > 3), the 
locations of zero amplitude can depend on both the locations and the strengths 
of the point sources, For a triple-point source, the first null points to occur 
as baseline length increases will satisfy the inequality. 



( 120 ) 


where 1 is either jPol or 1^,1, 
nr * 4 ' ' j ' 

To obtain fairly general results for structure effects for a triple point 
source, we will parameterize amplitude and phase as suggested by Eq. (IIA): 
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(1+g 


2 cos $2 + 


cos e^)' 


(g 2 sin 02 ■*’ ®3 


( 121 ) 


-g 2 sin ©2 ~ 83 sin 
tan iJijj l+g^ cos ®3 ®3 


( 122 ) 


where g. < 1. The BWS delay can be obtained as a function of the same four 
parameters by either numerically or analytically differentiating Eq. (122) with 
respect to w. For plots, the following numerical computation has been performed: 
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dui 


- .>3(83,63) 

em 


(123) 


where 

» (1 + e) 0 j^ 


(124) 
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and e is a very small number. An analytic expression for At, can also be 
readily 

? anted ^ . . 

> and f, in forming this BWS delay will be relative to the strongest point 
source. ’^To change the'^reference point for BWS delay to the brightness centroid, 
one must . et ?q in Eq. (106) equal to the centroid given by 


derived in terms of 83 * result will not be pre- 

here due to its uninformative c&inlexity. Due to the definition given to 


V 



i + $2 83 


(125) 


Ihe phase shift due to this change in Fq becomes 




(126) 


^ 2^2 

i ®2 ^3 


(127) 


Since (O^, O 3 ) are proportional to w, the BWS delay due to will be given by 



Ils. , i 82^2 ^ 3^3 

3 w w 1+82 '+• 83 


(128) 


Thus Che total structure delay becomes 




+ 



(129) 


wtem referenced to the brightness centroid, where At, and At^^ are obtained 
from Eqs. (123) and (128). 


Before proceeding to example plots, it is useful to derive an expression 


for the extrema in BWS delay as a function nf 82 ® 3 ' 

when + grt > 1 , the amplitude will go to zero for some values of (Qo»® 3 ^* 


As discussed above, 

and 

the BW§ delay will be infinite at those points. For + g„ < 1, the* 

BWS delay does not becotne infinite but peaks for certain values of ( 62 , 63 ). The 
extrenia in BWS delay at those points can be computed as follows. A consideration 
of phasors shows that, when g 2 + 83 < 1 , the most rapid change in as a func- 
tion of e„, 0 n occurs when • 
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(2m + l)ir 

m “ 0, + 1, +2, . . . 

(130) 

(2n + l)ir 

n » 0, + 1, +2, ... 

(131) 


The extrema of BWS delay will occur near these points. To obtain approximate 
values of those extrema, Eq. (129) can be analytically evaluated as a function 
of (Q2» 63) 3nd computed for the specific values of (62, 0^) In Eqs. (130) and 
(131). Such a calculation gives 


ATjj(m,n) 



g2(2m + 1) + g3(2n + 1) 
1 "* ^^2 ^ 3 ^ 


(132) 


where it is understood that §2 83 Figure 13 plots Axg as a function of 

(B2» 83) nt X-band for the first extremum, m = n = 0. From this example plot, 
one sees that Axg = 0.26 nsec for the particular case 81 “ g2 “ 'W'® 

saw for the case of a double-point source, structure effects can also be con- 
siderable for triple-point sources. 


For examples. Figs. 14 through 22 display amplitude, phase and BWS delay as 
a function of (0-, 0„) for two cases of (g,, g^)* 'the axes (p = 62/211 and 
q - 0 2/211) arc in units of cycles rather than radians and are equal to the pro- 
jections of the two separation vectors in the direction of the_^sky-projected 
baueline vector, with both in units of resolution (i.e,, B, • . Interpre- 

tation of these plots in terms of phasor sums such as those in *Tig. 11 will be 
left to the reader. To help clarify results, contour and surface plots of BWS 
delay are both shown. Even though the conversion of (0_, @2) (’^j v) might 

complicate geometric interpretation of these results, the values for the extrema 
in BWS delay show^ in_^the plots will not change with variable transformation or 
with changes in (P2» P2^ * Thus, given only point strengths (g2» g2)» one can 
readily determine for a triple-point source the extremal BWS delays (relative 
to the brightness centroid) either by inspecting plots such as Fig. 14 or by 
evaluating Eq. (132). 


As shown in Figs. 14 and 15, the structure delay for the case g2 = 0.2 and 
g„ = 0.4 passes through an extremum in delay whenever p and q are both multiples 
or 1/2, as explained earlier. The first extremum at p = q = 0,5 has a magnitude 
of about 0.1 nsec. A check on results for g2 -t g2 < 1 is supplied oy a compari- 
son of the extremal values predicted by Eqs. (130} - (132) with the results 
numerically computed to obtain these figures. The extrema in delay correspond 
to minima in fringe visibility shown in Fig. 16. These minima are all equal to 
about 0.25. Structure phase for this example is given in Fig. 17 and ranges 
between +0.1 cycles. 
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BWS DELAY RELATIVE TO BRIGHTNESS CENTROID, nsec 


10 



FLUX 2 

^2 " “flux 1 

Figure 13. Magnitude of First Relative Extremum in BWS Delay 
for a Triple-Point Source 
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Figure 15. Contour Plot of BHS Delay Relative to Brlghtneaa Centroid 
® for a Triple-Point Source, Case 1 
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Figure 18. BUS Delay Rel 


















The second example Is the case g„ " 0*4 and g« » 0.8. The delay values 
for this case are plotted In Figs, 18; 19 and 20, These figures are all plots 
of the same function and are -;iven together to help clarify the rather complicated 
behavior of delay. The surface plots in Figs. 18 and 19 are truncated above 1 
nsec to remove infinite spikes. When g 2 + Sq > tke amplitude will always 
be zero for some baseline vectors, as explained earlier. Null amplitudes In this 
example occur at adjacent pairs of points that are regularly placed in the p-q 
plane, as shown in Figure 21. Each pair corresponds to the + signs In Eqs. (117) 
and (118). For each null point, the delay approaches plus Infinity on one side 
of the point and minus infinity on the other, as shown in Figs. 19 and 20. Struc- 
ture phase for this example is shown in Fig, 22, The heavy jagged lines in that 
figure indicate an integer-cycle jump in phase from -0.5 cycle to 0,5 cycle. 

This discontinuous behavior has no ultimate significance and exists in this form 
due to the decision to keep phase between +0.5 cycle, as discussed earlier. 

Given that BWS delay is the partial of phase with respect to frequency, it is an 
interesting exercise to predict the plotted delay behavior on the basis of the 
structure phase plotted in Fig, 22, 


This example suggests that, if a source can be exactly represented by a 
triple-point source, it will be necessary to avoid certain critical values of 
baseline vector (provided the point sources are sufficiently far apart). Since 
the fringe visibility becomes very small at these points, it would be possible 
to alleviate this problem by placing a lower limit on allowed values of fringe 
visibility. 


With regard to observed sources, it is unrealistic to assume that a real 
source can be exactly modeled with three point sources. Actual distributions 
possess extended components rather than "point" components and/or have additional 
weaker components or background features. These characteristics make the truly 
singular behavior of discrete-point sources an unusual occurrence. Nevertheless, 
real sources can possess regions in which the fringe visibility becomes very 
small and the structure delay passes through extrema of considerable magnitude. 
Examples of real sources are given in the next section. 


XIV. EXAMPLES OF MEASURED BRIGHTNESS DISTRIBUTIONS 

The preceding sections have treated the analytical cases of a symmetric 
source, a double-point source and a triple-point source. Although analytical 
examples are instructive, only a thorough study of many actual source distributions 
can give a complete picture of source structure effects. To begin such an in- 
vestigation, this section reports on an analysis of ten brightness distributions 
measured by the GIT VLBI group. For each distribution, fringe visibility, fringe 
phase, BWS delay and effective position have been computed as a function of the 
sky-projected baseline vector (u, v), with the last three effects specified rela- 
tive to the ordinary centroid. Although instructive, this set of sources is not 
necessarily a typical set since sources specifically selected for structure 
investigations are often chosen on the basis of unusual characteristics. Another 
pctential difficulty is that the u-v coverage in these measurements is not as 
complete is some u-v regions as one would like. In such cases, computed structure 


effects must be used with caution* particulotXy BWS delays. 


The brightness disti itions resulting from interferometry measurements can 
be represented in a number of forms. As indicated in Appendix A* one representa- 
tion Is a set of delta functions placed on a regular grid. For this representation* 
the derivation for an N-point source in Section XIII (see Eqs 105-113) can be 
directly applied but the equations for phase* delay, amplitude and effective 
position should first be expressed as a function of the usual (u, v) coordinates. 
This is easily accomplished by rewriting Fq. (110) in the form 


\ - 2it(uBj^ + VYj^) 


(133) 


whore (6. , y. ) are the plane-of-the-sky coordinates of the^k point, expressed 
in radians, ^A similar expression can be written for the «Bg • 

(107), Further, the delay expression (Eq, 112) can be rewritten as 


where 



(134) 


(135) 

(136) 


The vector Pq is the reference position while ^P^^^ is the effective position of 
the source, where both positions are expressed relative to structure coordinates. 

In this work, the reference position Pq will be set equal to the ordinary centroid, 
as computed on the basis of the measured distribution. The quantity f is the 
observing frequency in hertz. Once a given source has been specified by a flux 
array t (g|<;, 3^, Y^^) , k ■ 1,N], the aforementioned equations can be used to compute 
phase, amplitude, delay and effective position. In the present computations of 
BWS delay and effective position, the frequency dependence of .the distributions 
will be neglected. 


Table 1 summarizes the 10 brightness distributions measured by the Caltech 
group (Refs. 2 and 3) and gives the source name, date of the measurement, observ- 
ing frequency, maximum values of u and v, total flux, interferometer stations, and 

a quantity k to be explained below. Because the total number of plots is excessive, 
s 
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Table 1. Summary of Brightness Distribution Heasuresaents 
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Station code: 1 — IJRAO; 2 — FDVS; 






complete results for only one of the measured distributions will be presented 
here. Brightness contours for the example distribution (9 in Table 1) are given 
in Fig. 23. The structure effects for this distribution are given in Figs. 24a-e 
in the form of contour plots as a function of u,v. In the plots, u is defined to 
be positive to the east and v positive to the north, Since the u-v coverage 
associated with a given distribution should not extend beyond the maximum allowed 
by baseline length considerations, each contour plot is marked with an approximate 
boundary outside of which the results are deleted. One Important point concerning 
this analysis is that the BWS delay has been computed on the basis of an ''artifi- 
cial" frequency of 8.3 GHz (see f in Eq. 134). In effect, this assignment of f 
pretends that, for BWS delay computation, the distribution was measured for the 
specified u-v values at 8,3 GHz (A ** 3.6 cm), even though it was not. The reason 
for this assignment is that most of our work will be carried out at X-band and it 
is therefore important to obtain BWS delay results at that frequency. To obtain 
the BWS delay at the actual frequency, one can easily scale the results by f“l 
(or by A) as in Eq. (134). All of the other plotted quantities require no expli- 
cit frequency assignment in computation, although they actually pertain, of course, 
to the observing wavelength given in Table 1 for distribution 9. 


The results for the 10 distributions were fairly complex but some general 
descriptive statements can be made. The magnitudes of the structure delays (the 
BWS delays computed for X-band) relative to the centroid were as large as a 
nanosecond but typically were less than 150 psec. The largest delays (*' 1 nsec) 
occurred in very localized regions in the u-v plane where very small fringe vlsi' 
bilities (~ 0.03) occurred. On average, structure delay increased as the fringe 
visibility decreased. 


With regard to effective position for BWS delay and phase-delay rate, the 
position of some sources changed as a function of u,v by as much as 10-30 marcsec 
while other sources remain fixed to within a few milllarcseconds. Thus, in 
applications that require position determination at 10 marcsec level, one must 
be careful in selecting sources and/or in correcting for source structure. For 
the specific case, of AVIBI applications, a final assessment cannot be made until 
intercontinental structure measurements (preferably at S-band and X-band) are 
available. As indicated in Table 1, the distributions analyzed here were produced 
by transcontinental baselines and therefore have inadequate resolution for inter- 
continental applications. 


The process of calibrating source positions for short-baseline measurements 
through the use of positions obtained on long baselines can also be investigated 
on the basis of the effective position plots. A more direct approach, however, 
is to study the structure delay plots along the anticipated u-v tracks for the 
proposed baselines. In general, the hounds on .structure delay discussed above 
also apply to this calibration process. 


XV. AN UPPER LIMIT FOR THE SOURCE STRUCTURE EFFECT IN BWS DELAY 

The preceding section has shown that source structure can sometimes contri- 
bute errors in BWS delays in the 10-30 cm range. Thus, if centimeter-level 
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Figure 23. Measured Brightness Distribution for 3C345, Distribution 9 
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Figure 24e. Effective Position (milliarcsec), Y Component, Relative to 
Brightness Centroid for 3C345, Bistributlon 9 




baseline measuremencs are a goal, some method must be devised for reducing or 
calibrating structure effects. A possible method is calibration through the use 
of measured brightness distributions, as discussed in the preceding sections. 

The primary difficulty with this approach is that an individual VLB! structure 
measurement is currently an expensive and time-consuming process. The prospect 
of working with very large catalogues, possible containing many tbno-varylng 
members, makes the calibration approach even less Inviting. It is therefore 
worthwhile to ins^estlgate alternate methods for overcoming structure problems. 


Another approach is suggested by the general tendency of BWS structure delay 
to increase with decreasing fringe visibility. Suppose a general formula could 
bo established that, purely on the basis of the value of fringe visibility, sets 

an upper limit on structure effects in BWS delay. Then, if the limit turned out 

to bo sufficiently small for some upper range of visibility values, the larger, 
unacceptable structure delays could be eliminated by merely deleting observations 
with the smaller visibility values. Such an approach is attractive since fringe 
visibility can often be obtained along with eacii VLBI observation. Thus, if the 
visibility determinations were accurate enough, the experiment would not depend 
on outside measurements for structure corrections. One important point that 
can be made is that the proposed upper limit would not have to be an absolute 
limit, valid for every source. For example, it would be useful to establish, 

if possible, on approximate 3cr statistical limit so that structure effects could 

be treated like other errors. Another example would be a limit that was valid 
for all sources except for Infrequent pathological cases. In fits with redundant 
observations, the pathological cases could be discovered and deleted through 
residual analysis. 


To begin development of a limit approach, this section derives and tests a 
limit formula for BWS delay due to structure. First, a fairly general analysis 
is used to derive an expression for an upper limit on structure delay. Next, 
special analytical sources are used to generate another expression for the upper 
limit. These theoretical results are compared with results from 10 measured 
distributions. 


The general derivation of a delay limit proceeds as follows. Since the 
results for a double-point source suggest that possible frequency dependence in 
a given brightness distribution will produce a fairly small de)'ay effect relative 
to that produced by the frequency dependence of the sky fringes, we will ignore 
the former dependence in the following derivation. As shown in Fig, 25, plot a 
given value of the brightness transform R in the complex plane along with a vector 
indicating the change in R produced by a change in frequency Aw. The quantities 
Z and Z are the real and imaginary components of R defined in Section III. The 
associated BWS structure delay will be given by 



Aw 


(137) 


where A4'g is the change in phase of the brightness transform. Among the possible 
directions for AR, the maximum value for Ax„ occurs when AR is perpendicular to R. 
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Figure 25. Schematic Illustration of the Change in Structure 
Phase with Frequency 


71 




Thus, an upper limit for Aig is given by 




(DZ^/Su)^ + (3Z^/3»)^ 

2 2 
2+2 
c s 


T 1/2 


(138) 


where the inequality holds if AR is perpendicular to R. This expression can be 
further reduced for a frequency-independent N-polnt source. For such a source, 
the frequency derivatives can be obtained from Eqs. (108) and (109), which give 


3w 




(139) 


(140) 


where we have used the fact that 8j^ is proportional to w. If we momentarily 
neglect w, the numerator in Eq. (138) becomes 


Numerator = 




sin 0j^ 


2 


1/2 


(141) 


This expression can be viewed as the length of a vector formed by summing N 
vectors of length 8i,|0}^l* ^’he length of the sum vector can be no greater than 
the sum of the lengths of the component vectors. Thus Eq. (138) becomes 


ATj < 



(142) 


For most real sources at larger u-v values, this limit will be much larger than 
the actual structure delay since the point-source vectors will add destructively 
to a large extent. 
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Up to thl»? point, the derivation has involved no approKlmations and placed 
no restrictions on the structure of the sources. Unfortunately, It now becomes 
necessary to make assumptions since it is possible to concoct a source that will 
generate any arbitrary set of 0. values. Assume a valid limit can be obtained 
by pulling outside of the sum a^factor 2Trk that represents an upper limit for 
the overall effective value of the 6j^ terml, This operation gives 


'‘■'b - t (J, 2 + z 2) 


which can be rewritten as 


CATjj 


5 k 


s V. 


(144) 


where we have expressed the results in dimensions of length and have used the 
definition of fringe visibility's 


V as 

s 


(z ^ + z 

c s ^ 


(145) 


This result indicated that the upper limit for structure delay varies in inverse 
proportion to fringe visibility and in direct proportion to observing wavelength. 
However, the limit is not of quantitative value unless an estimate for k^ is 
obtained . 


The results for measured brightness distributions in Section XIV can be used 
to make an estimate of k . The simplest procedure is to determine for each dis- 
tribution the value of k that makes the limit equation an absolute limit. This 
can be carried out by computing for each u-v point the value of k required to 
make the limit equation give the BWS delay obtained at that point from the bright' 
ness distribution. The resulting k^ values can be searched over the allowed 
region of the u-v plane to find a maximum value for' each distribution. The 
results are given in Table 1, which shows that the maximum values for k fell in 
the range of 0.11 to 0.31. 


Four interesting features became evident in the analysis of this particular 
limit equation. First, the u-v points at which the limit value was reached were 
usually not at fringe visibility minima although they were associated with Che 
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lower range of visibility values. Second, these points occurred near the u-v 
boundary for eight out of the 10 brightness distributions while even larger (but 
necessarily disallowed) values occurred outside tliat boundary. This result suggests 
that some of the largest allowed values (see Table 1) may be merely an artifact of 
"boundary effects." More information is required to test this possibility. Third, 
for the ijnportant visibility values in the range 0.1 - 0.2, as well as for larger 
visibility values, the actual BWS delays were almost all less than about half the 
values predicted by the limit equation with kg * 0.25. Thus, for these sources, 
the RMS delay error at a given visibility value was much less than the limit value, 
perhaps 1/3 as large. Fourth, for u-v points near zero, the delay effects were a 
very small fraction of the value predicted by Eq. 144. This result is not unexpect- 
ed since the BWS delay is actually zero at u * v » 0, while the limit formula gives 
k^X. This last faatura shBg'^sts that a better limit equation should be found. 


Another limit formula can be obtained from an apparently unrelated result for 
a double-point source. AC the first fringe visibility minimum, the BWS delay 
effect (see Eq. 94) can be written as 



(146) 


where the strict definition of Vo in this case is the value of the fringe visibility 
at the minimum,. Although Eq, (146) was not derived as a limit, we will test it as 
a limit by letting v be any given value of visibility. In comparing the two 
limits in Eqs. (144)°and (146), we find they are approximately the same at small 
visibilities (v^ s 0.2) if we set ■ 0.25. Since this is the approximate value 
for ks obtained above through comparisons with measured distributions, we see that 
Eq. (146) is approximately valid as a limit at small visibilities for those dis- 
tributions. As visibility approaches 1.0 near u » v ■ 0, Eq. (146) correctly 
approaches zero, unlike Eq. (144). As discussed in Appendix B, one can be even 
more thorough and show that Eq, (146) sets a valid (though loose) upper limit for 
a range of small u-v values. The encouraging results for these two regimes 
support the use of Eq. (146) as a possible limit formula. 


Although it is not conclusive, more support for Eq. (146) as a limit formula 
can be gained from theoretical considerations. First, the limit in Eq. (146) 
actually applies to a class of sources. For example, the result in Eq. (132) for 
a particular class of triple-point source can be reformulated at the first 
visibility minimum (m =* n * 0) so that it matches Eq, (146). Further, although 
is will not be shown here, a particular class of N-point sources can also conform 
to this equation. Second, despite the fact it is based on special cases, Eq. (146) 
provides a more general limit than is immediately obvious. For both of the limit 
steps (i.e,, Eqs. 138 and 142) established in the preceding derivation, these 
special cases require the equal sign. That is, AR is perpendicular to R for Eq. 
(138) and, for Eq. (142), all phasors resulting from the frequency derivative of 


-1 ^^k\ 

the point-source phasors (i.e., w coherently. Thus, the full vector 

strength of the derivative phasors is realized, as represented in the numerator of 
Eq, (142). Although there are these vector advantages, one can point out that the 
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derivative nraplltudos are limited since the 0^ values never exceed ir for the 
special cases. Thus for larger u-v values where the [Oj^j values are quite large, 
Che limit might fail. Whether failure occurs very often with real sources must 
be tested with many measured brightness distributions. One factor working against 
failure for many sources is the likelihood that increases in the 0j<; values will 
be counteracted by rather severe deviations from the coherent addition of phasor 
derivatives. 


At least two categories of distributions can be specified that would violate 
the limit formula. The first category consists of sources with a large, strong 
diffuse component displaced by a considerable distance from a weaker compact com- 
ponent. In this case, the effective position would move from the centroid of both 
components to the center of the compact component as baseline length increased. 

If the separation of components were great enough, large changes in structure delay 
could occur. The second category includes any source with two components that 
are separated by a distance much larger than the extent of either component. As 
baseline length increased, such a source could be modeled approximately as a 
double-point source until one of the components began to resolve out. Since many 
"sky fringes" could be placed between the components, the relative maxima for BWS 
delay in Eq. (94) might be present for several values of n. Since the limit form- 
ula is based on the first and smallest relative maximum (n * 0) , the larger maxima 
would exceed the limit formula, As mentioned earlier, the existence of such non- 
conforming sources docs not necessarily lnvalid£.<^;a the limit approach. As long 
as the limit formula only falls for a small percentage of sources, it still would 
be of great value in reducing and estimating structure effects in BWS delay for 
the other sources. As an example of the limit approach, suppose one set a lower 
limit of 0,2 on allowed fringe visibilities. Then the maximum structure delay 
predicted by the limit equation would be about 4.5 cm at X-band, as shown in Fig- 
ure 26. Hany observations would have visibilities larger than 0.2 and would 
therefore have smaller structure delay limits. Further, as discussed above, 
actual structure delays for the analyzed sources were almost always much less 
than the limit, with the RMS value being about 1/3 the limit at a given value of 
visibility. If we hypothesize on the basis of this preliminary statistical infor- 
mation that the maximum delay of 4,5 cm is an approximate estimate of the 3-o 
delay for all observations, then the l-o structure delay for all observations 
would be about 1.5 cm. This calculation suggests one might be able to reduce the 
1-0 structure delay to about 1.5 cm simply by deleting observations with fringe 
visibilities less than 0,2. For shorter regional baselines, one would be able to 
place a higher cutoff on allowed visibility values and therefore could reduce 
structure delays even more. For example, with visibilities greater than 0.5, the 
trial limit formula and the same statistical assumptions indicate a l-o structure 
delay of about 0.5 cm at X-band. 


As mentioned in Section XIV, in some applications the source positions 
obtained from long baseline measurements are used as a priori for short baseline 
measurements. One Important feature of such position calibrations is that any 
structure errors in delay encountered on long baselines will be reduced in effect 
by about the ratio of baseline lengths. This feature allows relaxation of the 
structure delay requirement that one must place on the long-baseline behavior of 
sources to be calibrated for use with regional baselines. For example, if l-o 
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structure delays equal to about 5 cm corrupt the source positions determined on a 
6000'»km baseline, those position errors would show up as about a 0.5-cra error 
(l-o) in the calibrated delays on a 600- km baseline. Since the preceding calcula- 
tion indicates that a l-o structure delay less than 5 cm at X-band on a long base- 
line should not be difficult to realize, position calibration errors due to long- 
baseline structure can probably be restricted to 0.5 cm or less for regional 
baselines. The other structure error in such calibrations, that due to variations 
in effective position on the regional baselines themselves, might be reduced to 
about 0.5 cm by a visibility cutoff, as discussed above. 


The preceding examples of the limit approach are based on the analysis of 
only a few sources and must therefore be treated with caution. To conclusively 
test the validity of Eq. (146) or any other limit formula, many more measured 
distributions must be analyzed using more sophisticated techniques. Improvements 
in the analysis would include computation of the following quantities for each 
distribution: (a) the RMS structure delay for u-v points with visibilities greater 
than a variable lower cutoff, (b) the number of u-v points with visibilities 
falling within separate visibility internals (e.g., with 0.1 spacing) and (c) the 
maximum and RMS structure delay for each of these bins. The results could be 
used to improve the limit equation by empirically determining a new multiplicative 
function of visibility to better account for the general dependence of the limit 
on visibility. Further, such an analysis would allow a delay value from the limit 
equation to be more reliably connected with a typical or RMS value for delay at a 
given visibility. 


XVI. SUMMARY AND CONCLUSIONS 


This report investigates structure effects in VLSI measurements and includes 
(a) a theoretical background, (b) an approach to structure corrections based on 
measured brightness distributions, (c) two analytical examples of sources, (d) 
the results of an analysis of measured brightness distributions, and (e) a limit 
formula for reducing and estimating BWS structure delay. In the theoretical 
derivation, the concept of effective position is formally introduced. It is shown 
that, for a given observation of any source, the source can be analytically and 
conceptually replaced by a point source located at the formal effective position. 
When BWS delay and delay rate are the observables, one unique effective position, 
valid for both observables, can be specified for each observation. For a given 
source, structure corrections are carried out by applying to each observable the 
correction that results from moving the effective position for that observation 
to an assigned reference position (e.g., the origin of structure coordinates). 

One insight offered by this approach is that the unknown origin of structure 
coordinates for an independently measured brightness distribution can be absolutely 
placed on the celestial sphere by means of standard VLBI measurements and multi- 
parameter fits. 


An analysis of the examples of a double-point source and a triple-point 
source shows that the largest and most troublesome corrections to BWS delay (>1 
nsec at X-band) and delay rate will occur in regions of the u-v plane near points 
of zero or nearly zero fringe visibility (<0.03). Near these points, structure 


77 


phase undergoes rapid variations as a function of (u,v) and the BWS delay 
correction can even (theoretically) approach infinity. In the eKample of a 
double-point source, the first line of zero amplitude in the u-v plane can be 
crossed only when the sky-projected baseline length satisfies Bg/X > 0.5/S, 
where S is the point-source separation. For a triple-point source, the first 
null point can be. reached only when Bg/X ^ 0.3/Sjt, where is the separation of 

point source 2 or 3 from the strongest point source. The analyses for double- 
point and triple-point sources also show that their effective positions (for BWS 
delay) can lie far outside the extent of the source near u-v points with very 
smal. fringe visibilities. In general, similar behavior in effective position 
will probably be observed in the u-v regions near local deep minima in amplitude 
for most sources possessing such minima. With regard to possible frequency 
dependence of brightness distributions, it is shown that a source consisting of 
two points with dissimilar spectral indices can require non-negligibie additional 
corrections to BWS delay at X-band (~150 psec « 4.5 cm) near points of very small 
fringe visibility (~0.03) if the difference in the two spectral indices is large 
(Ak " 1.0). However, this effect is about an order of magnitude smaller than the 
delay effect due to the frequency dependence of the sky fringes mentioned above. 
For larger, more common values of fringe visibility (>0.2), the delay effect 
will be less than 1 cm. This result suggests that it may be possible in general 
to neglect calibration of this particular structure effect whenever the goal for 
total delay error is greater than 1-2 cm. Since this hypothesis is based on 
results for a double-point source, the frequency dependence of many actual dis- 
tributions should be analyzed to see if the hypothesis is generally valid. Such 
a study has not been carried out to date. 


To begin a study of structure effects for real sources, 10 brightness distri- 
butions measured by the Caltech group have been analyzed. Again, the largest BWS 
delay effect (—1 nsec = 30 cm at X-band) occurred at points in the u-v plane where 
fringe visibility fell to very small values (~0.03). More typically, BWS struc- 
ture delay fell in the range of 1-10 cm. These results indicate that some method 
must be developed for handling source structure when delay accuracies of a few 
centimeters are desired. Effective position (for BWS delay) in some cases varied 
by as much as 10-30 marcsec, which indicates that, in applications concerned with 
0.01” position measurements, one must be careful in sei.ecting sources and/or 
correcting for structure effects. 


Since structure effects can be quite large, this report has begun an investi- 
gation of a limit approach to source structure effects in BWS delay. A trial 
limit formula (Eq. 146) has been constructed that might correctly set an upper 
limit on BWS structure delay. Theoretical considerations show that the formula 
sets a valid (though loose) upper limit for all sources as u and v approach zero. 
At present, the formula has been compared with too few real sources to establish 
its general validity. However, the seven real sources analyzed at this time sug- 
gest that the limit is approximately valid and that it is approximately equal to 
three times the RMS structure delay at a given visibility. If these results turn 
out to be valid for nearly all real sources over all u-v values, then the limit 
approach might be of great value in both reducing and estimating structure effects 
in BWS delay. For example, in the last section, rough calculations based on 
these results indicated that the 1-a structure delay at X-band might be reduced 
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uo about 1,5 cm by avoiding observations with visibilities less than 0.2 (e.g., 
on intercontinental baselines) and to about 0.5 cm by avoiding observations with 
visibilities less than 0.5 (e.g., on regional baselines). The structure delays 
calculated in these two examples,' are less than the total error budgets for liWS 
delay in some VLBI projects. Thus some projects might avoid, at least temporarily, 
the large potential burden Imposed by the present difficulty in obtaining VLBI 
brightness distributions. 


A final word of caution would be that the limit approach in its present form 
has not yet been verified. It may turn out that, at the large u-v values associ- 
ated with intercontixiontal baselines, too few natural sources possess sufficiently 
large fringe visibilities in the required u-v regions, F\»rther, the measured dis- 
tributions considered here are too few in number and too limited in u-v coverage 
and observing frequency to provide a general verification. The u-v coverage is 
too limited because (a) the absence of short baselines (< 300 km) can allow 
Important large-scale structure to be resolved out and missed and (b) the absence 
of intercontinental baselines can allow important small-scale structure to remain 
unresolved. Finally, in the present limit approach, there are two statistical 
assumptions that must still be generally verified: (a) that structure delays as 
an ensemble can be statistically modeled and (b) that the limit equation is equal 
to about three times the IIMS structure delay at a given visibility. In the last 
section, two categories of sources were specified that would violate the trial 
limit formula. More survey data is required to determine whether the percentage 
of sources that fall in those categories, or other disallowed categories, is too 
large for the limit approach in its present form to be valid. Even if the present 
approach is invalid it should be possible to construct another visibility-dependent 
limit that is valid for nearly all sources. How useful the final form will be in 
overcoming source structure problems remains to be determined. As a possible 
refinement to the limit approach, future studies should determine whether the 
nature of visibility signatures along a few u-v tracks can help eliminate unde- 
sirable sources. For example, one might show that sources with fairly smooth, 
constant signatures generally have smaller structure delays than those with widely 
varying signatures. To answer these questions, we plan to investigate additional 
sources in the near future and apply more detailed and sophisticated analysis, 
with the hope that the present limit approach, or some refined version, can be 
generally verified. 
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APPENDIX A. THE MEASUREMENT OF BRIGHTNESS DISTRIBUTIONS 


The determination of a brightness distribution for a given source by VLBI 
involves the measurement of fringe amplitude and fringe phase for that source 
for many baselines, with the various baselines preferably spanning a wide range 
of lengths. Measured amplitudes and phases (actually closure phase) are passed 
through an inversion procedure to determine the brightness distribution. For 
this inversion, the CIT group used the "hybrid-map" approach (Ref. A). 


In VLBI measurements, the spatial resolution can be no better than about 
X/B, where X is the observing wavelength and B is the maximum length among the 
baseline vectors. Consequently, without loss of Information, a measured distri- 
bution can be described by specifying a flux strength at each point on a grid in 
the plane of the sky where the dimensions of each square are set equal to a 
fraction of the best resolution. The output of the hybrid-map procedure specifies 
only the nonzero values and consists of a sequence of three numbers - a flux 
strength followed by two coordinates specifying the relevant grid point. Since 
the grid representation is unrealistically discontinuous and possesses more 
resolution than the Interferometer merits, the final hybrid map for display is 
computed by convolving the grid distribution with a Gaussian restoring beam whose 
dimensions are approximately equal to the measurement resolution. This operation 
smooths out the distribution so that the fine detail is more consistent with 
the resolution limits of the interferometer. Since it gives a more accurate 
visual Impression, the smoothed distribution will be shown for each source in 
this report. 


In the analysis to obtain structure phase, etc., one has the choice of 
working with the grid distribution or the smoothed distribution. An Important 
fact concerning this choice is the theorem that the Fourier transform of a con- 
volution is equal to the product of the Fourier transforms of the two convolved 
functions. Since the Fourier transform of the Gaussian restoring beam has zero 
phase, the structure phase values resulting from the two representations will be 
identical, as will the BWS delay and effective position. With regard to amplitude, 
the transform of the grid distribution will possess amplitudes that closely 
approximate the measured values along the u-v tracks of the interferometer. 

However, the transform of the smoothed distribution will possess amplitudes that 
are smaller than the measured values along the u-v tracks for the longest base- 
lines of the Interferometer. Since one would want the transform amplitudes to 
equal the measured values, the grid distribution will be analyzed in this report. 


Since sampling in the (u,v) plane is always incomplete, a measured bright- 
ness distribution is only approximately equal to the true distribution. Sampling 
in the (u,v) plane is sometimes incomplete in a local region. One example occurs 
near u = v = 0, which can cause large-scale structure to be missed. As mentioned 
above, sampling is always incomplete above the maximum (u,v) values, which means 
that the most compact components will not be evident. This report will not con- 
sider such deficiencies but will merely use the measured distributions to. mechan- 
ically compute structure effects for (u,v) points below the maximum values. 
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The problem of incomplete sampling in the u-v plane might be alleviated 
in some cases by the inclusion of four observables not currently used in structure 
measurements; BWS structure delay 0 <|>b/9w), phase rate 0<|)3/9t) and the fre- 
quency and time derivative of fringe visibility (9vg/3w, 3vg/3t). The phase 
derivatives would be separately summed around each baseline triangle to reveal 
structure information,- while the visibility derivatives could be obtained for 
each baseline. Measurement of the two derivatives with respect to frequency 
would require multichannel observations, such as those described in the text for 
bandwidth synthesis. 


To demonstrate how these new observables might complement the standard 
observables of visibility and phase, first consider the phase derivatives and 
rewrite BWS structure delay in the form 
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and phase-delay rate in the form 
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(A3) 



(A4) 


where 7 denotes the gradient with respect to u-v coordinates. In the equations, 
we have only considered the time and frequency dependence due to the sky fringes 
and have neglected the time and frequency dependence of the brightness distri- 
bution. As discussed in the text, results for a double-point source suggest that 
the variation of the sky fringes with frequency might produce an effect in the 
derivative that is about an order of magnitude greater than that due to frequency 
dependence of the brightness distribution (compare Eq. (91) with Eq. (94) and 
see p. 78), A study of a variety of real sources will be required to resolve 
this issue. On the other hand, the effect of the time dependence of the bright- 
ness distribution on phase rate is expected to be negligible. With regard to 
the utility of these observables, the phase rate observable might help to 
extrapolate into data gaps along a u-v track or might improve sampling in regions 
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of rapidly varying phase. The BWS delay observable might help to extrapolate 
into unsampled regions between tracks in the u-v plane since it is sensitive to 
the change in phase in the direction of such regions. The primary question is 
whether the accuracy of these closure observables, particularly BWS structure 
delay, will be sufficient to provide significant new information. Study of that 
issue Is beyond the scope of this report. 


With regard to the visibility derivatives, one can substitute Vg for in 
the above equations and show that the gradient of visibility in u-v coordinates 
could be determined by measuring 3<{)jj/Du and 3^ij/3t. We again assume that the 
time and frequency dependence due to the sky fringes greatly exceeds the time 
and frequency dependence of the brightness distribution. This is expected to be 
a valid assumption for the derivative with respect to time. Results for a 
frequency-dependent double-point source suggest that the frequency dependence of 
a source will contribute a much smaller effect to the visibility derivative than 
the sky-fringe effect, provided differences in spectral Indices of the components 
of the source are less than about 0.3. Again, more work is required to resolve 
this question for real sources. The extrapolation benefits described above for 
the phase derivatives also pertain to the visibility derivatives. To obtain 
accurate measurements of 3vg/Sm, one would have to make accurate relative cali- 
brations of the visibilities measured across the RF band. Again, it is beyond 
the scope of this report to estimate possible accuracies for the observables. 


One strong indication of the positive contribution that could be made by 
accurate determination of the derivatives of visibility and phase is provided by 
ordinary sampling theory. If one attempts to obtain a complete measure of a 
function on the basis of a uniform sample rate, then the measurement of both the 
value and the derivative of the function at each sample point will allow a re- 
duction of the minimum, saniple rate to 1/2 the sample rate required for value-only 
sampling. That is, the sample rate can be reduced from the Nyquist rate 
(Rfj = 2 X highest frequency) to Rjg/2. Thus, accurate measurements of the deriva- 
tives of visibility and phase would substantially Increase the effectiveness of 
a given set of inadequate u-v tracks if the contribution due to the frequency 
dependence of the distribution were sufficiently small for the source being 
measured . 
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APPENDIX B. FRINGE VISIBILITY, BWS DELAY AND PHASE FOR SMALL U-V VALUES 


For small u-v values, the 
visibility in Eq. (Ill) can be 


01^ values in Eq. (110) become small and the f ringe 
shown to approach 
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where 
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We have assumed that the 0^ values (which are expressed in radians) are evaluated 
relative to the ordinary centroid so that 
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Thus, for small u-v values, the fringe visibility decreases by an amount equal to 
the RMS spread of angles of the point-source phasors. Further, one can show that 
the BWS delay derived from Eqs. (112) and (113) can be reduced to the form 


3 

where < 9^ > is given by an expression similar to Eq. (B2) and where the reference 
position is again the ordinary centroid. 


The BWS delay result in Eq. (B4) can be used to test the limit formula in 
Eq. (146) for small u-v values. Although it was not derived for small 0^, Eq, 
(146) is used in limit applications by merely substituting a given value of 
visibility. Thus, that limit can be connected with the 0^ values through the 
use of Eq. (Bl). This operation gives 


cAt = I ^ 0^ > (^5) 

Due to the factor of it and due to the fact that the 0^ values are all assumed to 
be a small fraction of a radian, one can readily see that the actual delay given 
by Eq. (B4) will be much smaller than the limit predicted by Eq. 146 (i.e. , Eq. 
B5). Thus, at small u-v values, Eq. (146) is a valid upper limit for BWS delays. 
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since this limit is not as tight as one might desire, the form of Eq. (b 4) 
suggests a new limit formula that might be tested for small u-v values; 

b 

chx , Vg s 0.9 (B6) 

where the constant b is to be empirically determined but should be greater than 
1,0. This form will be tested along with other limit formulas in future work. 


Finally, a similar derivation can be carried out for phase in Eq, (107) to 

give 


A = -r < 0? .■> (in radians) (B7) 

P D K 

This result suggests a limit formula for phase delay that should ’^e tested for 
small u-v values; 

where b was defined above. Since the delay in Eq. (B8) is 1/3 as large as the 
delay in Eq. (B6) , this result for phase delay is another indication that phase 
delay can be less susceptible to structure than BWS delay. 
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